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Abstract. We consider a dynamics of Higgs' field in the framework of cosmological 
models involving a scalar field ^ from Nonsymmetric Kaluza-Klein (Jordan-Thiry) 
Theory. The field & plays here a role of a quintessence field. We consider phase 
transition in cosmological models of the second and of the first order due to evolution 
of Higgs' field. We developed inflationary models including calculation of an amount 
of inflation. We match some cosmological models, calculating a Hubble parameter and 
an age of the Universe. 



1. Introduction 

We develop in the paper cosmological consequences of the Nonsymmetric Kaluza- 
Klein (Jordan-Thiry) Theory. The theory has been extensively described in the first 
point of Ref. [1] . We refer to the book on nonsymmetric field theory for all the details 
of our theory and we present only cosmological applications of the theory. 

Let us give a short description of the theory. 

We develop a unification of the Nonsymmetric Gravitational Theory and gauge 
fields (Yang-Mills' fields) including spontaneous symmetry breaking and the Higgs' 
mechanism with scalar forces connected to the gravitational constant. The theory is 
geometric and unifies tensor-scalar gravity with massive gauge theory using a multidi- 
mensional manifold in a Jordan-Thiry manner. We use a nonsymmetric version of this 
theory. The general scheme is the following. We introduce the principal fibre bundle 
over the base V = E x G/Gq with the structural group H, where E is a space-time, G 
is a compact semisimple Lie group, Go is its compact subgroup and H is a semisimple 
compact group. The manifold M = G/Gq has an interpretation as a "vacuum states 
manifold" if G is broken to Go (classical vacuum states) . We define on the space-time 
E, the nonsymmetric tensor g a @ from N.G.T., which is equivalent to the existence of 
two geometrical objects 

g = g{ a p)T®e p (i.i) 

9=9[ a( 3]0 a Af (1.2) 

the symmetric tensor ~g and the 2-form g. Simultaneously we introduce on E two 
connections from N.G.T. W p~ and -Tg . On the homogeneous space M we define the 
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nonsymmetric metric tensor 




(1.3) 



where £ is the dimensionless constant, in a geometric way. Thus we really have the 
nonsymmetric metric tensor on or V = E x G/Gq. 



r is a parameter which characterizes the size of the manifold M = G/Gq. Now on the 
principal bundle P we define the connection to, which is the 1-form with values in the 
Lie algebra of H. 

After this we introduce the nonsymmetric metric on P right-invariant with respect 
to the action of the group H, introducing scalar field p in a Jordan-Thiry manner. 
The only difference is that now our base space has more dimensions than four. It is 
(ni + 4)-dimensional, where n\ = dim(M) = dim(G) — dim(Go). In other words, we 
combine the nonsymmetric tensor 7ab on V with the right-invariant nonsymmetric 
tensor on the group H using the connection lo and the scalar field p. We suppose that 
the factor p depends on a space-time point only. This condition can be abandoned 
and we consider a more general case where p = p(x,y), x £ E, y £ M resulting 
in a tower of massive scalar field pk,k = 1,2... . This is really the Jordan-Thiry 
theory in the nonsymmetric version but with (ni + 4)-dimensional "space-time" . After 
this we act in the classical manner. We introduce the linear connection which is 
compatible with this nonsymmetric metric. This connection is the multidimensional 

analogue of the connection on the space-time E. Simultaneously we introduce 
the second connection W. The connection W is the multidimensional analogue of the 
VF-connection from N.G.T. and Einstein's Unified Field Theory. Now we calculate the 
Moffat-Ricci curvature scalar R(W) for the connection W and we get the following 
result. R(W) is equal to the sum of the Moffat-Ricci curvature on the space-time E (the 
gravitational lagrangian in Moffat's theory of gravitation), plus (ni + 4)-dimensional 
lagrangian for the Yang-Mills' field from the Nonsymmetric Kaluza-Klein Theory plus 
the Moffat-Ricci curvature scalar on the homogeneous space G/Gq and the Moffat- 
Ricci curvature scalar on the group H plus the lagrangian for the scalar field p. The 
only difference is that our Yang-Mills' field is defined on (ni + 4)-dimensional "space- 
time" and the existence of the Moffat-Ricci curvature scalar of the connection on the 
homogeneous space G/Gq. All of these terms (including R(W)) are multiplied by some 
factors depending on the scalar field p. 

This lagrangian depends on the point of V = E x G/Gq i.e. on the point of the 
space-time E and on the point of M = G/G$. The curvature scalar on G/Gq also 
depends on the point of M. 

We now go to the group structure of our theory. We assume G invariance of the 
connection uo on the principal fibre bundle P, the so called Wang-condition. According 
to the Wang-theorem the connection oj decomposes into the connection cDe on the 




(1.4) 
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principal bundle Q over space-time E with structural group G and the multiplet of 
scalar fields Due to this decomposition the multidimensional Yang-Mills' lagrangian 
decomposes into: a 4-dimensional Yang-Mills' lagrangian with the gauge group G from 
the Nonsymmetric Kaluza-Klein Theory, plus a polynomial of 4th order with respect 
to the fields ^, plus a term which is quadratic with respect to the gauge derivative 
of <P (the gauge derivative with respect to the connection loe on a space-time E) plus 
a new term which is of 2nd order in the <P, and is linear with respect to the Yang- 
Mills' field strength. After this we perform the dimensional reduction procedure for 
the Moffat-Ricci scalar curvature on the manifold P. We average R(W) with respect 
to the homogeneous space M = G/Gq. In this way we get the lagrangian of our theory. 
It is the sum of the Moffat-Ricci curvature scalar on E (gravitational lagrangian) plus 
a Yang-Mills' lagrangian with gauge group G from the Nonsymmetric Kaluza-Klein 
Theory plus a kinetic term for the scalar field ^, plus a potential V(<P) which is of 4th 
order with respect to <P, plus £ int which describes a nonminimal interaction between 
the scalar field <P and the Yang-Mills' field, plus cosmological terms, plus lagrangian 
for scalar field p. All of these terms (including R(W)) are multiplied of course by some 
factors depending on the scalar field p. We redefine tensor g^ v and p and pass from 
scalar field p to \P 

p = e~*. (1.5) 

After this we get lagrangian which is the sum of gravitational lagrangian, Yang-Mills' 
lagrangian, Higgs' field lagrangian, interaction term £ mt and lagrangian for scalar field 
plus cosmological terms. These terms depend now on the scalar field <P\ In this 
way we have in our theory a multiplet of scalar fields (<P r , <P) . As in the Nonsymmetric- 
Nonabelian Kaluza-Klein Theory we get a polarization tensor of the Yang-Mills' field 
induced by the skewsymmetric part of the metric on the space-time and on the group G. 
We get an additional term in the Yang-Mills' lagrangian induced by the skewsymmetric 
part of the metric g a p. We get also Ci nt , which is absent in the dimensional reduction 
procedure known up to now. Simultaneously, our potential for the scalar-Higgs' field 
really differs from the analogous potential. Due to the skewsymmetric part of the 
metric on G/Gq and on H it has a more complicated structure. This structure offers 
two kinds of critical points for the minimum of this potential: $[? rt and ^rt- The first 
is known in the classical, symmetric dimensional reduction procedure and corresponds 
to the trivial Higgs' field ("pure gauge"). This is the "true" vacuum state of the theory. 
The second, ^l rt , corresponds to a more complex configuration. This is only a local 
(no absolute) minimum of V. It is a "false" vacuum. The Higgs' field is not a "pure" 
gauge here. In the first case the unbroken group is always Go- In the second case, it is 
in general different and strongly depends on the details of the theory: groups Go, G, 
H, tensors £ a b, g^ and the constants (, £. It results in a different spectrum of mass 
for intermediate bosons. However, the scale of the mass is the same and it is fixed 
by a constant r ("radius" of the manifold M = G/Gq). In the first case F(#° rt ) = 0, 
in the second case it is, in general, not zero V(<I>l lt ) / 0. Thus, in the first case, the 
cosmological constant is a sum of the scalar curvature on H and G/Gq, and in the 
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second case, we should add the value V(#J rt ). We proved that using the constant £ 
we are able in some cases to make the cosmological constant as small as we want (it 
is almost zero, maybe exactly zero, from the observational data point of view). Here 
we can perform the same procedure for the second term in the cosmological constant 
using the constant (. In the first case we are able to make the cosmological constant 
sufficiently small but this is not possible in general for the second case. 

The transition from "false" to "true" vacuum occurs as a second order phase 
transition. We discuss this transition in context of the first order phase transition in 
models of the Universe. The interesting point is that there exists an effective scale of 
masses, which depends on the scalar field 

Using Palatini variational principle we get an equation for fields in our theory. We 
find a gravitational equation from N.G.T. with Yang-Mills', Higgs' and scalar sources 
(for scalar field &) with cosmological terms. This gives us an interpretation of the 
scalar field W as an effective gravitational constant 

G eff = G^e-("+ 2 ^. (1.6) 

We get an equation for this scalar field \P. Simultaneously we get equations for Yang- 
Mills' and Higgs' field. We also discuss the change of the effective scale of mass, m e s 
with a relation to the change of the gravitational constant G e ff • 

In the "true" vacuum case we get that the scalar field \P is massive and has Yukawa- 
type behaviour. In this way the weak equivalence principle is satisfied. In the "false" 
vacuum case the situation is more complex. It seems that there are possible some scalar 
forces with infinite range. Thus the two worlds constructed over the "true" vacuum and 
the "false" vacuum seem to be completely different: with different unbroken groups, 
different mass spectrum for the broken gauge and Higgs' bosons, different cosmological 
constants and with different behaviour for the scalar field The last point means that 
in the "false" vacuum case the weak equivalence principle could be violated and the 
gravitational constant (Newton's constant) would increase in distance between bodies. 

We explore Einstein A-transformation for a connection W A g ((m+4)-dimensional) 
in order to find an interpretation of M + gauge invariance of W M field. We discuss M.+ and 
U(l)p invariance. We decide that U(l)p invariance from G.U.T. is a local invariance. 
Due to a geometrical construction we are able to identify W ^ from Moffat's theory 
of gravitation with the four-potential corresponding to the U(l)p group (internal 
rotations connected to fermion charge). In this way, the fermion number is conserved 
and plays the role of the second gravitational charge. Due to the Higgs' mechanism 

is massive and its strength, H^ v , is of short range with Yukawa- type behaviour. This 
has important consequences. The Lorentz-like force term (or Coriolis-like force term) in 
the equation of motion for a test particle is of short range with Yukawa-type behaviour. 
The range of this force is smaller than the range of the weak interactions. Thus it is 
negligible in the equation of motion for a test particle. We discuss the possibility of 
the cosmological origin of the mass of the scalar field \P and geodetic equations on 
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P. We consider an infinite tower of scalar fields coming from the expansion of 

the field \P(x,y) on the manifold M = G/Gq into harmonics of the Beltrami-Laplace 
operator. Due to Friedrichs' theory we can diagonalize an infinite matrix of masses 
for <Pfc transforming them into new fields &' k . The truncation procedure means here to 
take a zero mass mode $o and equal it to & . 

We discuss in the paper cosmological models involving field \P which plays a role 
of a quintessence field. We find inflationary models of the Universe and discuss a 
dynamics of the Higgs field. Higgs' field dynamics undergoes a second order phase 
transition which causes a phase transition in an evolution of the Universe. This ends 
an inflationary epoch and changes an evolution of the field \P. Afterwards we consider 
the field & as a quintessence field building some cosmological models with a quintessence 
and even with a -fT-essence. A dynamics of a Higgs field in several approximations gives 
us an amount of an inflation. We consider also a fluctuation spectrum of primordial 
fluctuations caused by a Higgs field and a speed up factor of an evolution. We speculate 
on a future of the Universe based on our simple model with a special behaviour of a 
quintessence. 

The paper is divided into three sections. In the second section we discuss infla- 
tionary models and phase transitions obtained due to dynamics of Higgs' fields. In the 
third section we present an evolution of Higgs field and quintessential models of the 
Universe. We consider several approaches to get an amount of inflation. Eventually 
we calculate an amount of an inflation in a simplified model of Higgs' field evolution 
and a power spectrum of primordial fluctuations caused by a Higgs field. 

2. Inflationary models and phase transitions 

According to new observational data [2] concerning distances of type la supernovae 
it seems that we need some kind of "dark energy" which drives the evolution of the 
Universe. This "dark energy" can be considered as a cosmological constant or more 
general as cosmological terms in field equations (in the lagrangian). In the case of cos- 
mological model this type of "dark energy" — "vacuum energy" is a cause to accelerate 
the evolution of the Universe (a scale factor R(t)) (see Ref. [3]). The cosmological 
constant is negligible on the level of the Solar System and on the level of the Galaxy. 
Moreover it can be important if we consider even nonrelativistic movement of galaxies 
in a cluster of galaxies (see Ref. [4]). In some papers considered cosmological terms 
result in changing with time of a cosmological constant (see Ref. [5]). Some of them 
introduce additional scalar field (or fields) in order to give a field-theoretical description 
of such an evolution of a cosmological "constant" . Those scalar fields are independent 
in general of the additional scalar fields in inflationary models. 

Thus the inflation field (or fields in multicomponent inflation, which can be 
the same as some of Higgs' fields from G.U.T. -models) can be different from those 
fields. In particular considering scalar-tensor theories of gravitation results in so called 
quintessence models (see Ref. [6]). Moreover in such theories there is a natural field- 
theoretical background for an inconstant "gravitational constant" . In such a way this 
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quintessence field can be used in twofold ways. First as a source of change in space and 
time of a gravitational constant. Secondly as a source of cosmological terms leading 
to the model of quintessence and a change in time of a cosmological "constant". In 
our theory we have a natural occurrence of these phenomena due to the scalar field ^ 
(or p). Let us consider the lagrangian of our theory paying a special attention to the 
part involving the scalar field •P', i.e.: 



(~ gauge 
e -(n+2)*£ YM(A) + _ _ £k . n( v ^ 

(n-2)-P e (n-2)U> „ \ 

-^ V <*>--SS=- £| -<*' A >) <"> 



r / e ^a 2 P 



* 1 . 



We put c = h = 1. 

Now we rewrite the lagrangian (2.1) in the following form. 



L = i?(W0 - 87rGjve- (n+2) *L 



matter 



- STrGiv-CscaK^) + e n * ( ^-a*12(r) + ^ 

V pi r / 

where in L m atter we include all the terms from Eq. (2.1) with Yang-Mills' fields, Higgs' 
fields, their interactions and coupling to the scalar field <P\ 
The effective gravitational constant is defined by 

G eff = G N e~W 



2<P 



P\ (2.1a) 



in such a way that the lagrangian of the Yang-Mills field in L ma tter is without any 
factor involving scalar field \P '. If the scalar field IF is constant (e.g. & = 0) we can 
redefine all the fields in such a way that we get ordinary (standard) lagrangians for 
these fields. 

Let us consider the situation after a spontaneous symmetry breaking and simplify 
to the case of g^ = 0. We get 

L = R - 87rG eff L matter - 8ttG n L sc ^) + SitG N U{V) (2.2) 

where 

£ SC ai = Mg^9 v ■ # >7 , M > (2.3) 
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8nG N U(V) = -Yl^i*^-*)* + e ("+2)* (^P^j + ^e"*. (2.4) 
We get the following equations 

~ ^ matter seal 

R»v ~ 2 R 9v = S^Ceff T^ v + 8itG n T^, (2.5) 

where and i? are a Ricci tensor and a scalar curvature for a Riemannian geometry 
generated by g^ v = g (piy) , 

matter 

T " F = ( P + ri//- P f (2.6) 
is an energy-momentum tensor for a matter considered as a radiation plus a dust. 

seal /M « \ — 

SvrGV 2^ = SvrGV f — ^ • ( 5 Q ^, a • «P i/3 ) - ■ $,„ J + <?^A C ^ (2.7) 



where 

777.— n- 2 ~ ~ ^n— ~ 

A c * = 2e^-All x V{^ rt ) - e^-^R(r) - e«*-Ap, (2.8) 



where mj is a scale of a mass for massive Yang-Mills' fields (after a spontaneous 

symmetry breaking for a true vacuum case). It is convenient to write 

_ e (™-2)<? e n<f e (n+2)<f 

A c ^ = - —7 ^— /3 (2-9) 

= P(S,a a ,m x ) (2.9a) 
7 = 7(C>n»2. a») 
K = 0, 1 corresponds to true and false vacuum case, i.e. 

F« t ) = 0, V(<Pl rt ) + 0, (2.10) 

a = 0, ai / 0. (2.10a) 
For a scalar field ^ we have the following equation: 



16irG N Mg al3 (v a (d^ + (n - 2)e^-^ a K - (n + 2)e ( - n+2 ^ 

- ne n * 1 - (n + 2)G eff T = 



(2.11) 
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where T = p — 3p is a trace of an energy-momentum tensor for a matter field. For we 
are interested in cosmological models we take for a metric tensor a Robertson- Walker 
metric: 

dr2 + r 2 d9 2 + r 2 sin 2 6d<p 2 



ds z = dt z - R\t) 



1 — kr 2 



, k = -1,0,1 (2.12) 



and we suppose that W, p,p are functions of t only. 
One gets 



V = JL h# - ( ^) a K e^* + i^±^/3 e ("+ 2 )* 
M 2 \ M J M 

"' ' (2.13) 



M 2 pl ^ 



+ 5 e " 7+ ^ e " (n+2) "- (/, - 3p) 



M p i - ' ■ — ' ' 



/8vr 

In this case we get standard equations for a cosmological model adapted to our 
theory 

^ = -8vrG eff Q(p + 3p)j - 8ttG n Q(^ + 3p*)) , (2.14) 

RR + 2R + 2k = 87rG eff Q(p~p)) # 2 + 87^ Q(p» -Pa*)) (2-15) 
Using Eqs (2.14-15) one easily gets 

ii + K = p-\ pq, (2.16) 



R 

where H = — is a Hubble constant 
R 



M ■ o 1- 
8ttG nP v = SnG N —y 2 + -A c k 

M -o 1- 
8irG N pq = 8irG N —V 2 - -X cK . 



(2.17) 



Let us consider a cosmological model for a "false vacuum" case, i.e. without matter 
and only with a scalar & and a vacuum energy V(<P\ rt ) / 0. In this case p = p = 
and we get 



6MR ■ , ~ dA c i 
# - M 2 , — — 

1 (l M 



2M<P —4r-M 2 1 -£- = (2.18) 



R 2 + k = - I _-^-^ + -A c i I R z (2.19) 
3 \ 2 M 2 j 2 i v 7 

^ = -Mj* + 2 A - (2 - 20) 



Let us take 
Thus we get 

and 



W = & ! = const. 
dXd 



(Mr 



■0?i) = o 



(n + 2)(3xj + njxl - (n - 2)«i = 
R 2 + k = ^X cl (x 1 )R 2 



Xi 



One gets from Eq. (2.25) 
where 



R(t) = Roe Hot 



A c i(xi) 
6 



is Hubble constant (really constant). From Eq. (2.23) we obtain 



2.21) 
2.22) 

2.23) 
2.24) 

2.25) 

(*) 

2.26) 

2.27) 

2.28) 

2.29) 
2.30) 

„r - ;\'n^r + K'H - LMi 1 - Lri, *(» + - ■ (2.31) 

2(n + 2)v3p 

Thus we get an exponential expansion of the Universe. Using Eq. (2.24) we get also 
k = (i.e. a flatness of a space). 

In this way we get de Sitter model of the Universe. This is of course a very special 
solution to the Eqs (2.18-20) with very special initial conditions 



(2.32) 



Xi 



-n7 + v /n 2 7 2 + 4(n 2 - A)a\f3 



2(n + 2)/3 



For oti > 0, (3 > we get 



n-2 



A c i(xi) 
H = 



X 



2(n + 2) 2 /3 

(n-2)/2 



y/n 2 ^ 2 + 4(n 2 - 4)qi/3 > n|'y| 

n7 2 - 7 V / n 2 7 2 + 4(n 2 - 4)«i/3 + 4«i/3(n + 2) 

/ 1 V2 

- 7V^ 2 7 2 + 4(n 2 -4)ai/3 + 4ai/?(n + 2) 



fl(0) 


— Ro 






§<°> 


= HoRo 


#(0) 


= #i 






*< 0) 


= 



Let us disturb the solution by a small perturbation and examine its stability. Let 

R = R + 5R, \SR\,\SR\ <^R. 
One gets in a linear approximation 



One gets 



d 2 \ 

2M<p + 6MH <p + M 2 1 -^(V 1 )i P = (2.34) 

^ = -~Mip 2 (2.35) 

+ A; = ^M0 2 (i! 2 + RSR). (2.36) 



,.r -- "" 'sin ( + (2.37) 



M 

where (we take for simplicity M = 



if 

4 

M < Mn = - 

3 



n(n + 2)7 V / n 2 7 2 + 4(n 2 - 4)«i/3 - 4(n + 2) 2 (n - l)«i/3 - n 2 (n + 2)7 2 
n7 2 + 4(n + 2)ai(3 - -f^n 2 -/ 2 + 4(n 2 - 4)qi/3 



(2.39) 



Thus we get an exponential decay of the solution ip, i.e. a damped oscillation 
around \& = \P\. It means the solution & = $i is stable against small perturbations of 
initial conditions for i.e. 

#(0) = #i + V2 sin5 (2.40) 
^(0) = 99 (-f^osinJ + cos^) . (2.41) 

Making and (5 sufficiently small we can achieve smallness of perturbations of initial 
conditions. The exponential decay of the solution to Eq. (2.34) can be satisfied also 
for some different conditions than Eq. (2.39) if we consider aperiodic case. However, 
we do not discuss it here. 
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Thus we get 

< Im^ 2 < a 2 e~ 3Hot (2.42) 

where a 2 is a constant. 

From Eq. (2.35) one gets 

< SR < a 2 e~ iHot , (2.43) 

which means that SR ~ and SR = O(SR(0)). It means SR = SR(0) + SR(0)t. 
Moreover from Eq. (2.36) we get that k = due to the fact that SR is a linear function 
of t. In this way we get 

R(t) ~ S(t) + 5/2(0) + «*(0)t = + ^^^ j e* ot . (2.44) 

This means that the small perturbations for initial conditions for i? result in a 
perturbation of i? i n such a way that i?o is perturbed by a quickly decaying function. 
Thus our de Sitter solution is stable under small perturbations and is an attractor 
for any small perturbed initial data. One can think that this evolution will continue 
forever. However, we should remember that we are in a "false" vacuum regime and 
this configuration of Higgs' fields is unstable. There is a stable configuration — a "true" 
vacuum case for which ax = (K = 0). For a = the considered de Sitter evolution 
cannot be continued. Thus we should take under consideration a second order phase 
transition in the configuration of Higgs' fields from metastable state to stable state, 
from "false" vacuum to "true" vacuum, from <P l crt to #° rt . In this case the Higgs fields 
play a role of an order parameter. Let us suppose that the scale time of this phase 

transition is small in comparison to — and let it take place locally. We suppose locally 

a conservation of a density of an energy. Thus 

seal seal matter 

Too = Too + T 00 . (2.45) 

We suppose also that the scalar field & will be close to the new equilibrium (new 
minimum for cosmological terms) and that the matter will consist of a radiation only: 

S ^oo = ^A c i(xi) (2.46a) 

seal matter j q 

T 00 + T 00 = -A c o(x )+8vrG J v^2. (2.46b) 

For xo = e* is a new equilibrium point we have: 

^(*o) = (2.47) 
A c0 = -e^ n+2 ^f3 - e^ 7 (2.48) 
^ = - e »*(( n + 2)/fe» + n7). (2.49) 
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One gets 



and supposing 



thus 



a *0 



x 



n\l\ 
(n + 2)/? 



A c0 (^o 



7 < 
2|7l 



2i _i_ i ri 
2 ^ 77, 2 



/?t( n + 2)f 



9+i 



(2.50) 

(2.51) 
(2.52) 



From Eqs (2.46ab) one gets 



1 



167tGjv I (n + 2) 2 /3 L 
2| 7 |f+ 1 



n 7 2 - 7^n 2 7 2 + 4(n 2 - 4)ai/3 + 2a 1/ 3(ra + 4) 



/3f(n + 2)f 



n+2 n + 2 

n 2 7 2 



(2.53) 



, \ n + 2 n + 2 

(n + 2) — /?— 



This gives us matching condition for a second order phase transition and simultaneously 
it is an initial condition for a new epoch of an evolution of the Universe plus a condition 
&(t r ) = 0, R(t r ) = R(t r ), where t r is a time for a phase transition to occur. Notice we 
have simply \P / tfi. Thus we have to do with a discontinuity for a field ^. Let us 
consider Hubble's constants for both phases of the Universe 



Hi 



6 

A c0 (x ) 



2 + 1 n 2 



6 3/?*(n + 2)i 



(2.54) 



Summing up we get 



Hi + H\ 

R(t r ) = R(t r ) 

&l(t r )=Mtr) 
tf 2 / il 2 



R e +H ° tr 
-- 



(2.55) 



Thus we see that the second order phase transition in the configuration of Higgs' 
fields results in the first order phase for an evolution of the Universe. We get disconti- 
nuity for Hubble constants and values of scalar field before and after phase transition. 

Let us calculate deceleration parameters before and after phase transition 



RR 



q = - 



R 2 



q = -l 



(2.56) 
(2.57) 
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before phase transition and 

q = -l (2.58) 

after phase transition. 

Let us come back to the Eqs (2.14), (2.16), (2.13). Let us rewrite Eq. (2.13) in 
the following form 

m+^9 + M*^=0. (2.59) 
R p (M 

(Remember we now have to do with a radiation for which the trace of an energy- 
momentum tensor is zero.) 

Let us suppose that ^ ps 0. Thus Eq. (2.59) simplifies 

_2_ M ^+^ = (2.60) 



Mpj <& 



and we get the first integral of motion 

-^pt^ 2 + A c0 = -5 = const. (2.61) 

2Af£ cU 2 v ' 

~M^ 2 = M^(S - 2A c0 ). (2.61a) 
One gets from Eqs (2.14) and (2.16) 



3R 1 p r 1 ( 2M f , 



R ' e ( n + 2 )^ 2 I M 2 { 



^ - A c0 (2.62) 



Using Eq (2.61a) we get: 



3R 1 p r - 5 , . 



One can derive from Eqs (2.64-65) 



^-(RR) = -\5R 2 + l\ c0 R 2 - k (2.66) 
at 6 3 



A c o = -e^(/3e 2 ^+7). (2.67) 
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Let us take k = and 5 = and let us change independent and dependent variables 
in (2.66) using (2.67) and (2.61a). One gets 



W - 1)0 + y ((n + 4), 2 - (n + l)y - 2) | = ^(y 2 - l)/(y) 



where 
and 

It is easy to see that y 2 > 1 for 

~dt 



± 



M pl | 7 |^ 



vV(y 2 - 1). 



(2.68) 

(2.69) 
(2.70) 

(2.71) 



Moreover according to our assumptions ^ « and this can be achieved only if < 
y — 1 < e, where e is sufficiently small. 

Thus for further investigations we take z = y — 1, y = z + 1. One gets 



d 2 f rlf 
2z(z + 2){z + l) 2 ^ + (z + 1) ((n + 4)z 2 + (n + 7)z - (n + 3)) J 



AM 

— Z ( z + 2 )f(z)=0. 



(2.72) 



Let us consider Eqs (2.64-65) supposing k = and (5 = 0. 

After eliminating R (via differentiation of Eq. (2.65) with respect to time t) and 
changing independent and dependent variables one gets: 



d_ 

dy 



8vrG 



N- 



Pr 



1 



+ ^"(y 2 - 1) ) / 



,n+2 



where 



It is convenient for further investigations to consider a parameter 

4M 



r = 



In such a way we get 



d 2 f df 
~ l) d^ + V ^ n + A)lj2 ~ {n + 1)y ~ 2 ^~ " l)f{v) 



(2.73) 
(2.74) 

(2.75) 

(2.76) 
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and 



d 2 f 

2z(z + 2)(z + 1 f^2+( z + !) (i n + 4 ) z2 + ( n + 7 ) z " ( n + 3 )) 



dz 



(2.77) 



-rz(z + 2)f(z) = 0. 
One can transform Eq. (2.73) into 



d_ 

dy 



FAT- 



y n+2 1 2' 

Changing the independent variable from y to z = y — 1 one gets: 

Pr 3 



+ -</V - 1) / 



y™" 1 [(n + 2)y 2 -n] / 2 



(2.78) 



8ttGn 



■(z + l) n z(z + 2)\f 



(z + l) n + 2 2 
(z + 1)"" 1 [(n + 2)(z + l) 2 -n] / 2 



(2.79) 



Let us come back to Eq. (2.71) in order to find time-dependence of <P '. One gets 

d& 



or 



/ 
/ 



t-U 

■ = ±^Jr M P 1 

y/p e ( n + 2 )* - Me"* v^f 
da; t — ti 



x\fj3x n+2 — \^\x r - 



M, 



'M 



P i- 



(2.80) 
(2.80a) 



If n = 21, where I is a natural number, one gets for f3 > 

2p+l 







'^£©(-2)^ 

,fc=l p=0 



2p 



^x + 



i 1 ^ v^-ItI-v^v^x + H) ' 

2^3 n Vv / ^"l7l + V3(v^x + l7l), 



(2.81) 



± (*-*i) 



M p i. 



In this case it can be expressed in terms of elementary functions. For a general 
case of n, (3 and 7 one gets 



2x-»iF a (^-5,(l-5);^) (t _ tl) 



M 



(2.82) 



where ii*2(«) &j c; z) is a hypergeometric function. For a small z (around zero) one finds 
the following solution to Eq. (2.77) 



f(z) = C lZ 



(n + 7)(5n + 23) ; 

1 grvTi) — ^+ C2 l 14 



(n-1) 



(2.83) 
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Ci, C2 = const. Thus 



R(y) 



d(y-l)^- 1 



+ C 2 1 + 



(n + 7)(5n + 23) 



8(n+ 1) 



(y-i; 



(n-1) 



(y-i) 



(2.84) 



where y > 1 (but only a little). 

Let us make some simplifications of the formulae taking under consideration that 
y — 1 is very small. 

One gets 



f(y) = c 1 + 



(n-l) 



(y-i) 



c > 



i?(y) = i2i (^l + ^-^-l 
After some calculations one gets 



(2.85) 
(2.86) 



2(n + 5) 4 10r(ra + 2) 3 
(n + 4)(n-l) y + (n- l)(n + 4) y 
(8r - (n - l)(n + 4)) , 2r(n + 3) (r - n + 1) 

+ , w . - - y - -7-^ — + 



2(n- l)(n + 4) 



(n 2 - 1) a (n-l) 



(2.87) 



(taking into account that (y — 1) is very small). Let us make some simplifications in 
the formula (2.71) for y close to 1. We find: 



dy 



'2 



Vy- 1 " /3f v 7 ^ 



^-(t-tlJMp! 



and finally 



^(t-ti) 2 + i 



(2.88) 



(2.89) 



2M/3- 

such that y = 1 for t = t\. 

Let us pass to the minimum value for p,- with respect to y close to 1. Thus we are 
looking for a minimum of a function 



1 



8vrG 



-y 



2(n+l) 



2(n + 5) 4 10r(n + 2) , 

ry + 7 — ttt — —^y 



n [ (n + 4)(n-l) a (n-l)(n + 4)' 

+ (8f- (n- l)(n + 4))^ 2 2r(n + 3)^ + (r - n + 1) 



2(n- l)(n + 4) 



(n 2 - 1) y (n-l) 



87tGtv 



(2.90) 
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for such a y that is close to 1. 

Let us write y = 1 + 77, where r\ is very small and develop W±(l + rf) up to the 
second order in r/. One finds 

W 4 (l + ?7) - V 2 {rf) (2-91) 



and 
where 



V 2 (v) = arj 2 + br, + c (2.92) 
(-n 2 + 3n(20r - 9) + 4(32r - 31)) 



2(n + 4)(n - 1) 
(-n 3 - 4n 2 (7r - 3) + 5n(18r - 11) + 22(3r - 2)) 
(n + 4)(n 2 - 1) 



(2.93) 
(2.94) 



_ (-3n 3 + 2n 2 (9r - 8) + n(50r - 29) + 4(7f - 4)) 
°" 2(n + 4)(n 2 -l) (2 ' 95) 

Let us find minimum of ^2(77) for rj > such that V^r?) > 0. For a < 0, 6 < 0, 
c > we get 

V2(r/ min ) = (2.96) 

and 

,. in = ^. (2.97) 

In this case 

(1 + r/ min ) 2(n+1) ^4(l + r?min) > (2.98) 

and will be close to the real minimum of the function (2.90). Simultaneously 

f?min < 1- 

For 

Pr = (2-99) 



(where a is a Stefan-Boltzmann constant) we get T m ; n and we call it a Td (a decoupling 
temperature-decoupling of matter and radiation). Thus yd = 1 + rj m in- Ud is reached 
at a time td (a decoupling time) 

t d = t 1 + ^0-M ph (2.100) 

From that time (to) the evolution of the Universe will be driven by a matter (with 
good approximation a dust matter) and the scalar field The radius of the Universe 
at t = td is equal to 

1 

R{td) = Ri(± + j^-^vL?) 2 • (2.101) 
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The full field equations (generalized Friedmann equations) are as follows: 



2 



[r) = -r-»» + z{2jif + 2 A -J < 2 - 102 > 

3R_8irG eS 1 / M - 2 \ 

If " " 2 y " Ac °J (2 - 103) 

2M - 6M R- d\ c0 . , B , ni . 

" jjjr • ^ " ^ + (» + 2)8vr Pm G eff = 0. (2.104) 

Before we pass to those equations we answer the question what is a mass of the 

scalar field ^ during the de Sitter phases and the radiation era. One gets: 

-m\ = ^fi^(^i) (2.105) 

1 2M d& 2 y ' K J 

-ml = ^^{%). (2.106) 

2M d<P 2 y ' K ' 

The second question we answer is an evolution of the Universe during a period from 
t r to ti, i.e. when y is changing from y = ^j-^pi to y = 1. 
This will be simply the de Sitter evolution 

R(t) = Re Hlt (2.107) 

where 



Hl = ] J^pl. (2.108) 

This will be unstable evolution and will end at t = ti, i.e. for iP ~ and value of the 
field & will be changed to the value corresponding to y = 1, i.e. 

^ l) = IlnA. (2.109) 

Thus we have to do with two de Sitter phases of the evolution with two different Hubble 
constants Hi and H . In the third phase (a radiation era) the radius of the Universe 
is given by the formula 

R(t) = R^HoU+H^-U) f l+ _\l\ n+4 y _ ti)4 \ 2 / 2 no) 

W V 3M(n-l)/3" V ' ) y J 

and such an evolution ends for t = td, 

R{t d ) = R exp{H t r + H 1 (t 1 -t r ))-(l + j^^rg ijl ) . (2.111) 
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In both de Sitter phases the equation of state for the matter described by the scalar 
field ^ is the same: 

p# = -p# m (2.112) 

In the radiation era we get 

pq, = \pq,. (2.113) 

Now we go to the matter dominated Universe and we change the notation intro- 
ducing a scalar field Q such that 

Q = ^TT- ( 2 - 114 ) 
M p i 



In this case we have 



PQ = \Q 2 +U{Q) (2.115a) 
Pq = \Q 2 ~ U{Q) (2.115b) 



where 



HQ) = \x c o (^) = -f + l| 7 | e «3 (2.116) 

a = Mpl , 6 = M pl . (2.117) 

VM VM 

Let us write Eqs (2.102-104) in terms of Q: 

Q + 3HQ + U'{Q) - 8ir(n + 2)G eS p m = (2.118) 



where 



G eS = G N e (2.119) 

H 2 = i (svrGefrp™ + ^Q 2 + E/) (2.120) 

^ = 1 (87rG effPm + 2{Q 2 - uj) . (2.121) 
Let us define an equation of state 

Wq = — = \&~ U (2.122) 

Q PQ \Q 2 + U 1 ' 
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and a variable 

_ 1 + W Q _\q 2 

XQ - Y^W Q ~ —> (2 ' 123) 

which is a ratio of a kinetic energy to a potential energy density for Q. It is interesting 
to combine (2.118) and (2.120) using (2.122) and (2.123). 
Let us define 

-Q 2 + U 

^ = (2.124) 

n m = ^^. (2.125) 

One gets 

fi m + Qq = 1 (2.126) 

or 

87rG eSPm = 3(1 - Q Q )H 2 . (2.127) 

We call p m an energy density of a background and for such a matter the equation of 
state is 

W m = V — = 0. (2.128) 

Pm 

It is natural to define 

p m = 8irG eS p m (2.129) 

as an effective energy density of a background with the same equation of a state 
as (2.128). 

Under an assumption of slow roll for Q, i.e. 

l -Q 2 « U(Q) (2.130) 

one gets 

W Q ~ -1, (2.131) 

i.e. an equation of state of a cosmological constant evolving in time. In this case 

xq ~ 0. This can give in principle an account for an acceleration of an evolution of 
the Universe if we suppose those conditions for our contemporary epoch. 
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3. Evolution of Higgs Field and Quintessential-Cosmological Models 

Let us come back to the inflationary era in our model. For our Higgs' field is 
multicomponent, i.e. it is a multiplet of Higgs' fields, we have to do with so called 
multicomponent inflation (see Ref. [7]). In this case we can define slow-roll parameter 
for our model in equations for Higgs' fields. 

One gets from the first point of Ref. [1] (see Eq. (5.7.7), p. 385) 



| (Lf)_ - 3H (Lf) m = -^S^L g A (3.1) 

v ft ) at) 



where L$~ = Lf~ (a time component of L^~) is defined by 

ldcL d Q ~ a + l cd garhg™L$ d = 2l cd g~ afh g™ c j t , (3.2) 
SV ( SV 1 

V, — F and <^ — r \ are defined in [1] (Eq. (5.6.9), (5.6.10), (5.424)). 

h y n J av 

3H 2 = 8nG N I ^ V(<P) - — — £ kin (<2>) + ^Ki^i) J (3.3) 

/ Op" 2 *! . \ 

H = 8nG N ^£ kin (<Z>) (3.4) 



where if is a Hubble constant and 

£ kin (#) = u{/"L« 5 ^|} . (3.5) 

<Zq is a constant value for a field & '. The slow-roll parameters are defined by 

>. > av 



m (3-6) 

, ^i( L u)i<A)L (37) 



with the assumptions of the slow roll 

e = 0( V ), 5 = 0( V ) (3.8) 

for some small parameter 77. 
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Under slow-roll conditions we have 



\ h / av 

and with a standard extra assumption 

A = Otf). (3.11) 

We remind that the scalar field & is the more important agent to drive the inflation. 
However, the full amount of time the inflation takes place depends on an evolution of 
Higgs' fields (under slow-roll approximation). The evolution must be slow and starting 
for <P = <Pl rt . It ends at <P = <P® It . 



Thus we have according to Eq. (2.26) 



^ ~ m UU 6nd \ ~~ H (t end — initial) 
-K (^initial J 



where 



<2>(t end ) = <2> c ° rt (3-12) 

^initial) = $lt (3-13) 
^end 

(we omit indices for (p). 

Thus we should solve Eq. (3.9) for initial condition (3.12) and with an assumption 
H = H . We have of course a short period of an inflation with H = H\ from i en d = t r 
to t = ti, i.e. 

JV tot = N + JVi = H (t r - t initial ) + Fx(t! - t r ). (3.14) 

Let us remind to the reader that N is called an amount of inflation. Let us consider 
inflationary fluctuations of Higgs' fields in our theory. We are ignoring gravitational 
backreaction for the cosmological evolution is driven by the scalar field $ '. We write 
the Higgs field as a sum 

^(f,t) = ^(t)+S^(r,t) (3.15) 

where ?^(t) is a solution of Eq. (3.9) with boundary condition (3.12-13) and <5#^(r, t) 
is a small fluctuation which will be written in terms of Fourier modes 

S0Ur,t) = J2^'kJ t y^- ( 3 - 16 ) 

K 
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The full field equation for Higgs' field in the de Sitter background linearized for ^"gAf) 
can be written 

d 2 f , r j \ „ „ d ( d \ 1 ?2 / / 



«• K Md ^^ ZH °f t { M h)^ir> K \ M i 



av 



#V'(&Atj) I ( 3 - 17 ) 



where 



l dc M%~ h + lcd9~a m g^Mi. = 2l cd g~ a7h g fnZ 5^ c . (3.18) 
Let us change dependent and independent variables. 

dr==£-, R(t) = Roe Hot (3.19) 
K(t) 

*M— -U (3.20) 
— oo < r < (a conformal time) (3-21) 

and 

a 

5$% . = -Ml . (3.22) 

Km ft v ' 

Simultaneously we neglect the term with the second derivative of the potential. Even- 
tually we get 



<i 2 



dT 2 

where 

^ ' i^i, : 'cti:/am:j ''/,';- -' ciiuinii .7 \ ^- 



l dc M%~ h + l cd gam9™Mi = 2l cd g~ am g™ c X % ? . (3.24) 



Let us notice that \Kt\ is the ratio of the proper wavenumber \K\R 1 to the 

Hubble radius — . At early times \Kt\ S> 1, the wavelength is small in comparison to 
Ho 

Hubble radius and the mode oscillates as in Minkowski space. However, if r goes to 
zero, \Kt\ goes to zero too. The wavelength of the mode is stretched far beyond the 
Hubble radius. It means the mode freezes. Thus the analyzes of the modes are similar 
to those in classical symmetric inflationary perturbation theory neglecting the fact that 
we have to do with multicomponent inflation and that the structure of representation 
space of Higgs' fields should be taken into account. 

However, the form of Eq. (3.23) is exactly the same as in one component symmetric 
theory if we take Mi„ a field to be considered. The relation (3.24) between Mi„ and 
X K ~ is linear and under some simple conditions unambigous. 
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It is possible to find an exact solution to Eq. (3.23) and Eq. (3.24). One gets 

a _ =r a (T\K\cOs(T\K\)-Sm(T\K\)\ 
^Km lKrh 



and 



+ C 2Km 



M d -~ = Ar~ 

Kb 1Kb 



~ 2Kb 





T 




t\K 


sin (ri?) 






r 




t\K 


cos (tK[ 






r 




t\K\ 


sin(r \K \) 


+ cos(r \K \) 



(*) 



(**) 



where 

In this way one obtains 



l dc N?~ + Ugamg^N?^ = 2l cd g~ aih g™Cf-„, i = 1, 2. 



-H n t 



«5<f>« = J_ e --ffo* y a ( _6 



and using (*) one easily gets 

\K\ 



Krh lKrh \ RqHq 

\K\ 



e~ Hot COS 



\K\ 



■e~ Hot + sin 



\RoH 



+ ^2/?™ I — I R(jHo 



K\ 

sin I ' „ e 



Hn * + cos 



\K\ 



RqH 



-H t 



(***) 



(****) 



(V*) 



where 



Ca tt /~ia 

■ i> - — — -"0^ ■ $■ - ■ 

i Km i Km 



(VI*) 

However, it is not necessary to use a full exact solution (V*) to proceed an analysis 
which we gave before. 

Thus the primordial value of Rg(t) is equal to 



R K = - 



(3.25) 



\t=t* 



where R K (t) is defined in cosmological perturbation theory as 



(3.26) 
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(see Ref. [8]) and t* is such that \K\ = R(t*)H Q , i.e. 



t* 



In 



\K\ 



RqHq 



\K\ = K. 



(3.27) 



The primordial value is of course time- independent. We suppose that fluctuations of 
Higgs' fields (the initial values of them) are independent and Gaussian. 

Thus we can repeat some classical results concerning a power spectrum of primor- 
dial perturbations. One gets 

2 



\ U Krh\ 



Ho 
2vr 



Thus 



Pr(K) 




(3.28) 



(3.29) 



t=t* 



for t* given by Eq. (3.27). 

Taking a specific situation with concrete groups H, G, Gq and constants £, £, r we 
can in principle calculate exactly P R (K) and a power index n s . Thus using some 
specific software packages (i.e. CMBFAST code) we can obtain theoretical curves of 
CMB (Cosmic Microwave Background) anisotropy including polarization effects. 

Let us give the following remark. The condition for slow-roll evolution for <P can be 
expressed in terms of the potential V. If those are impossible to be satisfied for some 
models (depending on G, G ,H, G' and parameters £, £), we can employ a scenario with 
a tunnel effect from "false" to "true" vacuum and bubbles coalescence. In the last case 
the time of an inflation (in the first de Sitter phase) depends on the characteristic time 
of the coalescence. 

Let us consider a more general model of the Universe filled with ordinary (dust) 
matter, a radiation and a quintessence. 
From Bianchi identity we have: 



d . .3R 

faPtot + (Ptot +Ptot)-^- = 



(3.30) 



where 



Ptot = PQ +Pm+Pr (3.31) 
Ptot=PQ+Pr- (3.32) 

Let us suppose that the evolution of radiation, ordinary matter (barionic + cold dark 
matter) and a quintessence are independent. In this way we get independent conser- 
vation laws 



(3.33) 
(3.34) 
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and 



One gets 



PQ + (PQ+P Q )^§ = 0. (3.35) 



^4 

p m = , -4 = const. (3.36) 

p r = — , B = const. (3.37) 
R 



From Eq. (3.35) we get 



p Q + (l + W Q )p Q -^ = 0. (3.38) 

Substituting Eqs (3.36-37) into Eqs (2.14) and (2.16), remembering Eqs (2.114) and 
(2.116-117) one gets 

{ il\ 2 , 1 (A B\ 1 



D +k= wiA : « + w) + u^r R2 (3 ' 39) 



IT " " (aw + aw j " 2MJ (1 + 3Wq)pq (3 ' 40) 

and from Eq. (2.118) 



Q +3H ®- { i^P = - (3 - 41) 

However, now Q and J7 (Q) are redefined in such a way that factor appears be- 
fore pq. This is a simple rescaling. We have of course 

p m = e- (n+2) *> m = e~ aQ Pm (3.42) 

p r = e - {n+2) *p r = e~ aQ p r . (3.43) 



Thus we get 



4 e aQ (3.44) 



Pr = ^e aQ . (3.45) 

Let us consider Eqs (3.39-41) and Eq. (3.35). Combining these equations we get 
the following formula 

(n + 2) A 



(l + W Q )p Q = 0. (3.46) 
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The one way to satisfy this equation is to put Wq = — 1. It is easy to see that if A = 
then Eq. (3.46) is satisfied trivially. The condition A = does not imply B = and 

VM 

we can still have 5/0. For such a solution Q = Q = const, and Qo = &o found 



earlier 



Thus 



M p i 



*** = = \ j n ' 7 ' = - f^) J n ' £ l ~ ■ (3.47) 
• / (n + 2)/3 a s \m p J ]j ( n + 2)R(r) 



-aQ -(n+2)<f _ n n+2 



or 



J \ n J \ IP 



e -«Q. = WfV!^WB£>y . (3.49) 



.W V n J \ \P 
Let us put Wq = -1 into Eqs (3.39-40). One gets 

1 (A B 



3 * W (3.51) 



i? Af2, V2i? 3 3i? 4 



From Eq. (3.50) one gets 



/ go" 



R*-kR 2 + ^L- R + 



B 

pi pi 

if A = 0, 



A (3.52) 



±i?f/7? =dt. (3.53) 



PQ E>4 J. E>2 i g. 



3MJ " " 1 3M 

Taking x = R 2 one gets 



/ q ^ = = ±2(t - 1 ) (3.54) 

\/3M^ X ~ kx + 3Mj 

and finally 
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where 



B 
— V- 

PQ 



For a flat case k = we find: 

± r * .V!» (t _ t0) . 

7 a/2/ + 1 ^ 
The last formula can be easily integrated and one gets 



(3.56) 



(3.57) 



R 



B 



(3.58) 



where we take a sign + in the integral on the left hand side of Eq. (3.57). 
Let us calculate a Hubble parameter (constant) for our model. One gets 



H 



R 
R 



\/3M 



pi 



B 



B 



(3.59) 



For a deceleration parameter one obtains 

RR (2sh 2 (it) - ch(u)) 



R 2 



where 



ch(u) sh 2 (it) 



B 



■(t-to). 



(3.60) 



(3.61) 



Thus we get a spatially flat model of a Universe dominated by a radiation which is 
expanding and an expansion is accelerating. 

Let us take k = and B = in Eq. (3.52). One gets 



±RdR 



PQ 
3W, 



= dt 



R 4 + 



3M 



-R 



pi 




R=?l— X X = R? 

PQ V A 




One finally gets 



x dx 



yjx{x* + 1) 3M pl 



pQ{t-t ). 



(3.62) 
(3.63) 

(3.64) 
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Let us notice that 



H z = 
RR 



R 
R 



3M p 2 



A jl \ ^ PQ 

r^ +p V > m; 1 



(2 Pq R 4 - A)R 
2{A + p Q Ri) 



(3.65) 
(3.66) 



Thus the model of the Universe is expanding and accelerating. 

Let us consider the integral on the left hand side of Eq. (3.64). One gets after 
some tedious algebra: 



/ 



^/x(x 3 + 1) 
598 



VV3 + 2 
12 



(^9 - 2>/3) In W + ^ \jl + V3 (6y/S + ll) iT 
^92^/4v¥+3 (l41\/3 + 272) + - 15) V598\J 16 + 9y/sj F 



(3.67) 



where 



W = 



_ V3 (16 - 7^3) / 2x + ^3 - 1 \ 
2 2 Ux-v^-lJ 

V 



+ 2\/6 - yfl 



A 



/ (2, + V3-l) V + 7 _ 4V 3 
^(2 X - y/S-l)J 



(2x + yg - 1) \ _ 
(2x - -v/3 - 1) / 2 

/ 



"(2g + V3-1)' 
(2x - y/3 - 1) 



II = II \ arccos 



(3.68) 



F = F arccos 



(3.69) 
(3.70) 



where Fl is an elliptic integral of the third kind and F is an elliptic integral of the first 
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kind 



u 

Hn,k)= f g . 2 (3.71) 



sin 99 



n(u,n,k)= [ dy == (3.73) 

J (1 — nsin 99) v 1 — k 2 sin 2 99 

n = V3 - 1 (3.74) 



and with the same modulus as -F(w, fe) 

A; 



.. _ 5 - 2^3 
" 13 



However, the most important condition for an existence of a physical solution is coming 
from the first part of a sum in the right hand side of Eq. (3.67). 
We should have 

W > 0. (3.75) 
In order to analyze condition (3.75) let us write W in the following way 



-20-^ 



+ ^^y + 2 V ^73 ] J(y + 7-4V3) (y - 1 + f ) 



= (3.76) 

y- 1 

where 

Obviously y > 0. 

We have two possibilities: 

I y > 1 (3.78) 

and 



-20 - ^fi + (16 - 7V5) y + 2V6 - v^y (y + 7 - 4^3) (y - 1 + ^) >0 (3.79) 
II y < 1 (3.80) 
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and 



_20 - ^ + ^-^ y + 2^6 - Vs^(y + 7 - 4^3) (y - 1 + 4) < °- ( 3 - 81 ) 
The first possibility can easily be solved: 



or for the second possibility 



II 



2-V3 



We have finally 



0.7916 < x < 



v / 3+ 1 



For 



we get 



y > 4.0612 (3.82) 
< y < 1. (3.83) 

(3.84) 
(3.85) 

x > (3.86) 

x < 3.072. (3.87) 
Let us notice the following: the function W has a finite limit at x = . Thus 



3^-5 1 

^ < z < ~- 

2 2 



v^+l 



we can remove a singularity at x = v ^ +1 . In this way we get 



0.7916 < x < 3.072 
3\/3- 5 



or 



< x < -. 
2 2 

This simply means that the solution exists only for 



or 



A I A 

0.7916, 3 /— < R < 3.072 ? — 

PQ V PQ 



2 V PQ 2 W p Q 



(3.88) 
(3.89) 

(3.90) 
(3.91) 
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Let us calculate a density of a quintessence energy: 



PQ = M 2 pl U(Qo) = -le^>((3e^-\ 1 \) 



n \ z m 



(3.92) 



n + 2) \a 2 s {n+1) J \m p ij \R(r) 
An energy density of a matter is equal to 



R 3 

and 



Pm = 4l eaQ ° ~ PQ eaQ ° ( 3 - 93 ) 



Tl + 2 



.... . e * _ m (^ " ^ W . (3.94) 



Let us consider the behaviour of the "effective" gravitational constant in our model 

n+2 



GW = G N e~^+ 2 ^ = G N e~ a ^ = G N (J-^)' 



Thus G e ff is a constant. For we are living in that model we should rescale the constant 
and consider G e s(Qo) (Eq. (3.95)) a Newton constant. Thus 



„+2 n+2 • . H+2 



(3.96) 



where Go is a different constant responsible for a strength of gravitational interactions 
in earlier epochs of an evolution of the Universe. It means we should write 

a M = | GNa: ^ (^) ■** * + | . . W (3 , 7) 

The obtained solution (i.e. (3.67)) is for IUq = —1. However, this is not an attractor 
of the dynamical equations. This is similar to the tracker solution of Steinhardt (see 
Ref. [9]). Moreover we cannot apply his method because our equation for a scalar 
field & (or Q) is different. It contains an additional term with a trace of the energy- 
momentum tensor (matter + radiation). Only with a radiation filled Universe our 
equation is the same (if we identify p r with Steinhardt density of a matter). In this 
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case we could apply Steinhardt tracker solution method and apply his criterion for a 
potential U{Q). In our case U(Q) is given by Eq. (2.116) 

U'(Q) = — ^ ((n + 2)(3e aQ - n| 7 |e &Q ) (3.98) 
2VM 

M 2 

U"(Q) = |L ((n + 2) 2 /?e aQ - n 2 | 7 |e b °) . (3.99) 

4VM 

The Steinhardt criterion consists in finding 

U"U _ ({n + 2) 2 (3e aQ -n 2 | 7 |e b Q) (f3e aQ - | 7 |e 6 °) |()() 
" (U') 2 ~ ((n + 2)(3e a Q -n\j\e b Q" 2 ' 

To get a tracker solution we need 

r>i. (3.101) 

In our case 

4/ 3| 7 | e (a+b)Q 

r = 1 — ,r 3.102 

((n + 2)/3e a( 3 -n| 7 |e b( 3) 

r . i - - 4 N = i _ -M -*-%Q = 1 _ 4| 7 | w 

(n + 2)/9 /3(n + 2) /3(n + 2) V ; 

Thus asymptotically we are close to r = 1. Moreover the Steinhardt criterion is 
only a sufficient condition. Our radiation filled Universe seems to be unstable due to 
appearing of a matter (a dust matter). Moreover a quintessential period of a Universe 
model has a severe restrictions of a value of a radius. Thus the solution with Wq = — 1 
and a matter independently evolving cannot evolve forever. 

In order to answer a question what is a further evolution of the Universe we come 
back to Einstein equations with a scalar field 'F and matter sources. Supposing as usual 
a Robertson- Walker metric and a spatial flatness of the metric we write once again a 
Bianchi identity 

T W 1 . . matter 

4t^ + -TPTZ~ {n+ ) T =0 (3.104) 

^ M 2 M 2 f 

where 

matter 

T ^ = p m u^u v (3.105) 
is an energy-momentum tensor for a dust matter. One gets 

d\ c0 ■ M ... M - 2 3R 

-| —iLr tZn& \T/ 2 

<W M 2 M 2 R 



(3.106) 
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We make the following ansatz concerning an evolution of a matter density 

Pm = A>e (n+2),? (3.107) 
po = const. (3.108) 

In that moment the scalar field & and the matter p m interact nontrivially. 

Using Eqs (3.107-108), (3.106) and an equation for a scalar field & one obtains 

P ° (n + 2) Po + ^=0. (3.109) 



Ml Ml ™ M* R 



Using Eq. (3.109) and Eq. (2.102) one gets 



d* _ V2M pl (n + l\ 2 3p 3 

r dx = ±^ ( t - to) (3.1H) 

7 xy/6 + 3/?x™+ 2 + 3 7 x™ 



where 



x = e^. (3.113) 
Using Eq. (3.111) and Eq. (2.102) one gets 

UnR=^±^ (3.114) 

dt M pl v 7 

or 

i? = J R e (S ^ )( *"*° ) . (3.114a) 

Let us consider Eq. (3.111) and let us suppose that x < 1 (and practically small). In 
this case we can simplify 

/ ^ - / , ^ (3.H5) 
J x v / 6 + 3(3x n + 2 + 3 7 x™ 7 xV* + 3 7 x™ 

and finally we get 

^r = —In (J-) - -J lLM pl (t -t ). (3.116) 

2n \3\-y\) n\ 3M p V ; V 7 
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Thus our prediction that x < 1 and is small has been justified. Using Eq. (3.116) one 
easily gets 



PV = -Z-7T + „ fK eX P 



3n 2 2x/3 
3n 2 2^3 



exp 



2* „ , 
3M pV 7 



2f „ . 
3M PV ; 



(3.117) 
(3.118) 



It is easy to see that 



P* *- 



(3.119) 



Thus in this model we have to do asymptotically (practically very quickly) with a stiff 
matter equation of state: 

p = p. (3.120) 

Simultaneously an energy density for a scalar field H/ is dominated by a kinetic 
energy which is practically constant and 



2^ " 3n 2 



(3.121) 



In this sense a dark matter generated by \P in the model is in some sense if-essence 
(not a quintessence). 

Let us calculate an energy density for a matter 



Pm = Po{ 77777 ) exp | \l -=M p i(i - t ) 



3|7l 



n V 3M 



(3.122) 



The effective gravitational constant reads 



3|7l 



(3.123) 



It is easy to write R as a function of 1^: 



n + 1 / — 



dx 



Xo 



xy/5 + 3(3x n+2 + 3-yx r 



(3.124) 



where 



x 



3|7l 



2n 



(3.125) 
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Eq. (3.124) can be easily reduced to the Eq. (3.114) using Eq. (3.111). 

Thus an energy density of matter goes to zero in an exponential way. The effective 
gravitational constant is growing exponentially. What does it mean for the future of 
the Universe? First of all the Universe will be very diluted after some time of such an 
evolution. Secondly, a relative strength of gravitational interactions will be stronger. 
This means that if we take substrat particles as cluster of galaxies then in a quite short 
time any cluster will be a lonely island in the Universe without any communication 
with other clusters. All the clusters will be beyond a horizon. On the level of a single 
cluster the strength of gravitational interactions will be stronger and eventually they 
collapse to form a black hole. In an intermediate time a gravitational dynamics on the 
level of galaxies' clusters will be governed by a Newtonian dynamics (nonrelativistic) 
with a gravitational potential corrected by a cosmological constant. The cosmological 
constant induces a positive pressure (a stiff matter) and will be important up to a 
moment of sufficiently large gravitational constant. After a sufficiently long time only 
a classical Newtonian gravity will drive a dynamics of galaxies' clusters. Let us roughly 
estimate this effect. In order to do this let us suppose that 

5M 2 

p& = p& = 2 = const. (3.126) 



In an energy momentum we get 



seal <5M 2 <5M 2 , 

T v = 2 ■ —^-u^u" - —^-g^. 



3n 2 



3n 2 



(3.127) 



However, we have in the place of Newtonian constant Gat a G e s which is a function 
of cosmological time and in the equation for gravitational field for g^ v we should put 



(M p f) 2 



seal 

rj-i {lis 



M, 



pi 



3ra 2 1 M p f 



(2«"u" - </H 



(3.128) 



where 



n+2 



eff 



pl 



m 



2(n + l) 2 - 3p 
3IPI 



x exp 



(3.129) 



Thus 



seal 



(M;ff 



T" v = Ae- K(t ~ to) {2u»u v - g» v ) 



(3.130) 
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where 



5 ym pl m~J y 3|£| J 

Consider the following case (the justification will be given below): 

'tV = -Ae- K{t ~ to) g^ . (3.133) 



(M p f) 

Now let us solve the following problem. Let A be a value Ae~ K ^~ to ^ in an established 
cosmological time t = t and the same for G = G e s(t). Let us consider a gravitational 
field of a point mass Mq static and spherically symmetric. Consider a metric 

ds 2 = B{r)dt 2 - -^j- - r 2 (d9 2 + sin 2 9 d6 2 ) (3.134) 

in spherical coordinates and Einstein equations 

R v - 2 9 ^ R = ~ A 9»v (3.135) 

for the metric (3.134) where 

A = Ae~ K{i - to \ (3.136) 
Via a standard procedure one obtains 

2GM Ar 2 

B(r) = l-— —-— ■ (3.137) 
Thus an effective nonrelativistic Newton-like gravitational potential has a shape 

V (r)= ««+^. (3 ,38, 

Remembering that G is growing exponentially in time and A is decaying in time ex- 
ponentially we see that the second term is important only for a short time. For a 
contemporary time 

Z~10- 52 m- 2 (3.139) 

according to the Perlmutter data and it could be important on the level of a size of 
galaxies' cluster. However, in a short time it will be negligible even on that level. 
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Let us calculate a Schwarzschild radius for a mass M with an effective gravitational 
constant G e g. One gets 



„eff 



2M G = J2M G eH (i) = r g 



a s J 



n+2 



n + 2 



3IPI 



2(n + l) 2 - 3p 



exp 



(3.140) 



Thus r| (calculated for an effective gravitational constant G e g(t)) is growing 
exponentially. It means that after some time it will be of an order of size of galaxies' 
cluster. Then galaxies' clusters collapse to form black holes. The time of a collapse is 
easy to calculate. 



Tg '"cluster 
Collapse = t + (n + 2)V^M 



(3.141) 



In 



^cluster 



(n + 2)ln( mp1 ^ 



3IPI 



2(n + l) 2 - 3p 



VS(n + 2)M pl ] (3.142) 



where r g is a gravitational radius of a mass of a galaxies' cluster with G = Gjy. 

Let us calculate a Hubble constant and a deceleration parameter for our model. 
One gets 



H = | = (n + l)M pl 
- - i 



(3.143) 
(3.144) 



Summing up we get the following scenario of an evolution. The Universe starts in a 
"false" vacuum state and its stable evolution is driven by a "cosmological constant" 
formed from P(P), P, F(#J rt ) and !Pi being a minimal solution to an effective self- 
interacting potential for <P This evolution is stable against small perturbation for 'F 
and P. The evolution is exponential and a Hubble constant is calculable in terms of 
our theory (the nonsymmetric Jordan-Thiry theory) . This evolution ends at a moment 
the "false" vacuum state changes into a "true" vacuum state. In that moment an en- 
ergy of a vacuum is released into radiation (and a matter). The Universe is reheating 
and a big-bang scenario begins at a very hot stage. Moreover the evolution is still 
governed by a scalar field \P which attains a new minimum of an effective potential. 
The effective potential is different and a new value of <p a is different too. The 
evolution is exponential and a new Hubble constant Hi is also calculable in terms of 
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the underlying theory. In this background we have an evolution of a multiplet of scalar 
fields from (p\ n to ^rt- This goes to a spectrum of fluctuations calculable in the 
theory. After that time the field & is slowly changing & ~ and the temperature of 
the Universe is going down. The evolution of a radiation is in the second inflation- 
ary stage adiabatic and afterwards during next phase (W ~ 0) nonadiabatic. The last 
means there is an energy exchange between radiation and a scalar field \P '. The total 
amount of an inflation N to t is calculable in the theory. The evolution of a factor R(t) 
is governed by a simple elementary function of a small rise. After this the radiation 
condenses into a matter (a dust with zero pressure) and a matter (a dust), a radiation 
and a scalar field evolve adiabatically without interaction among them. The evolution 
of a scalar field is governed by a quintessence, i.e. pq = WqPq, where Wq = —1. 

The quintessence field Q is a normalized scalar field During a radiation era 
~ 0) an equation of state for a scalar field ^ is different: pq, = (similar 
to stiff matter equation of state). The scale factor R(t) evolves according to some 
elementary function built from hyperbolic function in radiation + quintessence scenario 
and according to complicated function formed from the first and the second elliptic 
integrals (after reversing of those functions). It expands and accelerates. The solution 
with radiation and quintessence has a restricted behaviour (Eq. (3.90-91)). The energy 
density of radiation is going to zero. Only a quintessence energy density is constant. 
The energy of a quintessence is a potential energy dominated. However, the solution 
cannot be continued up to the end for an interaction between a matter and a scalar 
field starts to be important. The effective gravitational constant remains constant 
during the period. Let us remind to the reader that during inflationary epochs the 
effective gravitational constant was constant. Only during a radiation era it was slightly 
changing. In order to keep the Universe to evolve forever it is necessary to find a 
solution with an interaction between a matter and a scalar field $ '. Using an appropriate 
ansatz for an evolution of an energy density of a matter we get such a solution. The 
scale factor R(t) is exponentially expanding. 

An energy density of a matter is going to zero and an energy density of a scalar 
field is approaching (exponentially) a constant. However, now a scalar field & forms 
a different form of a matter — a stiff matter with an equation of a state p<p ~ p& 
(Wq, ~ 1). The effective gravitational constant is growing exponentially in time. In 
contradistinction with the previous period of an evolution the scalar field \P forms a 
matter with a kinetic energy dominance (i.e. -fT-essence) . In the previous period we 
have to do with a potential energy dominance (i.e. quintessence). 

Let us give some remarks on a spatial curvature of a model of the Universe. In the 
first period (inflationary scenario) the spatial curvature has to be zero and any pertur- 
bations of initial conditions cannot change it. In next periods it is natural to suppose 
that it is not changed. Thus our model of a Universe is spatially flat. In the moment of 
change of a phase of an evolution we should apply a matching conditions for an energy 
density and for a scale factor (a radius of a model of a Universe). However, an evolu- 
tion in such moments undergoes phase transitions for a law of an evolution of a scalar 
field and a scale factor, which can be considered as a second order or even a first order 
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phase transition. It is worth to notice that in formulae concerning Hubble constants, 
deceleration parameters, gravitational constants (effective) meet some parameters from 
unification theory of fundamental interactions and from cosmology. It seems that we 
are on a right track to give an account for large number hypothesis by P. Dirac. This 
demands of course some investigations especially in developing nonsymmetric Kaluza- 
Klein (Jordan-Thiry) theory. It is also interesting to mention that the scalar field \P 
plays many roles in our theory. It works as an inflation field during an inflation epoch. 
(However Higgs' fields play an important role in creating fluctuations spectrum.) It 
plays also a role of a quintessence and a iT-essence, except its influence on an effective 
gravitational constant and on an effective scale of masses in unification of fundamental 
interactions. This field is massive getting masses from several mechanisms. During 
inflation phases it acquires mass due to spontaneous symmetry breaking, different in 
both epochs. It gets a mass from a cosmological background. Thus except its cosmo- 
logical importance the fluctuations of a scalar field around its cosmological value could 
be detected (in principle) as massive scalar particles of a large mass. In this way an 
effective gravitational potential in nonrelativistic limit takes a form 

v{r) = ^ _ £ e _ ( ^ eK(t _ to) + ^! e - K(t - to) (3 _ 145) 

where k, A are given by Eqs (3.131-132), r is a range of massive scalar 5^ (a fluctuation 
of a scalar field & in a cosmological background). This could be detected as a tiny 
change of a Newton constant in gravitational law (e.g. the fifth force) on short distances. 
Let us notice that on short distances and on short time scales an effective gravitational 
potential reads: 

(l _£«-<*>). (3.146) 

On large distances and short time scales 

= G~Mo + *! 
r 6 



On short distances and large time scales 

G N M / _ a 
r \ r 

On intermediate distances and large time scales 



v{r) = ^ _ £ e -(^ gK(t - H (3148) 



V(r) = ^° e «(*-*o). (3 . 149) 

The latest being a cause to collapse a cluster of galaxies into a black hole. It seems 
also possible to consider an effective coupling constant of for a s enters and m e j?. 
However 



mf = e 



* (^) , (3.150) 
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r is a radius of a manifold G/Gq (e.g. a radius of a sphere S 2 ). If we rewrite Eq. (3.150) 
as 

pff 

-£ = £ e"*, (3.151) 

it would be quite difficult to distinguish a drift of rrf^ from a| ff drift. Only a quotient 
has a definite dependence on the scalar field W '. Especially it is interesting to consider a 
drift of Q em (a fine structure constant) reported from several sources. However, without 
extending our theory to include fermion (spinor) fields this is impossible. Thus it is 
necessary to construct a nonsymmetric-geometric version with supersymmetry and 
supergravity including noncommutative (anticommutative) coordinates to settle the 
problem. The scalar field & can form an infinite tower of massive scalar fields which 
could have important astrophysical and high-energy consequences. 

Let us consider Eqs (5.2.51) and (5.2.56) from the first point of Ref. [1], p. 318-319, 
and changing a base in the Lie algebra f) we get 

A = % (3.152) 

and 

a\ = 5 C % . (3.153) 
Simultaneously for G C H (q C h) we get 

fk = C k ( 3 - 154a ) 

fk = c k ( 3 - i54b ) 

where a, 6, c refer to the complement m in q (q = 0o+m), i, J, k to the subalgebra go 
of q. Let us suppose a symmetry requirement 



From Eq. (3.155) it follows that 



[m,m] C flo- (3.155) 



^ = 4 = 0- (3.156) 



Eq. (2.5.51) looks like 
Eq. (2.5.55) 



n°r« = ^QV (3-157) 

Let us come back to Eq. (3.1). 

For our inflation driving agent is a scalar field & we do not carry any backreaction 
of Higgs' fields and we can consider some nonstandard scenarios of an evolution of 
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those fields. Thus we can consider an evolution which is not a slow-roll evolution. Let 
us suppose that 

(3.159) 



— ($ a ~ ) ~ 0. 
dt { m> 



In this way we get from (3.159) 



J t ( Ld h)av + 



2r 2 \ 5$l 



9~bf, 



(3.160) 



The next step in simplification is such that we collapse all the degrees of freedom of 
Higgs' fields into one 

(3.161) 



<51 = Jjl d> 



In this way constraints (3.157) are satisfied trivially and Eq. (3.158) is as follows: 



H2r = Cir ($ 2 ~ 1 - 0) 

ab ab \ ) 

In this way V = V (no constraints) and 

A 



where 



and 



A 



l ab (2g^g^C^C b di - C^g™ g^L^ 



L a ~ = L a r ( a s <P 2 - — 
ab ab y as 

hc9rhig cm Lfa + l C dgamg mc Lf 5 = 2l cd g arh g mc Cf~. 
We can get the following identity 

l dc g^g^Ll d C c M + l cd g^g^l\fl~ q = 2l cd g^g^C(jO c M . 
Using Eq. (3.160) we get the following equation for <£: 

d 2 <P e™* 1 SV 



dt 2 + 2r 2 d$ 
if the following constraints are satisfied 



= 



(3.162) 

(3.163) 
(3.164) 

(3.165) 
(3.166) 

(3.167) 
(3.168) 
(3.169) 
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Thus from Eq. (3.168) one gets a first integral of motion 



<P 2 e"*" 1 B 
— + ——z-V = — - = const. 
2 2r 2 2 



From Eq. (3.170) one gets 



/ . - * -*(«-*,) 



r 2 a 2 



A{a 2 <P 2 -a s <P-iy + B 



and finally 

/ ± = ±V^(' - «.) 

•Po 

where 



a = a, 



5r 2 



= — y>. 

In order to find an integral on the left hand side of Eq. (3.172) and its 
function we use a uniformization theory of algebraic curves. Let 

f(<p) = (a-<p 2 + v + l)(a + f 2 - if - 1) = -<p 4 + 2(p 3 + <f 2 - lip + (a 2 - 1). 

One easily notices that 

1 + V5 + 4a 
^o = j 

is a root of the polynomial f(<p) (a real root). 
In this way one gets 

v -HP(z)-hf"M) +(P0 

where 

P(z) =P(z;g 2 ,g 3 ) 
is a P Weierstrass function with invariants: 

5 -3a 2 

52 = g 

a 2 10 

53 " 12 + 27 • 
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One easily finds 



f'(<Po) = |(5 + 3a)V5 + ^ (3.181) 
f"(<po) = -2(5 + 6a) (3.182) 



and 



In this way one gets 



z 

'<Pa 



r[f(x)]-Ux. (3.183) 

J inn 



W 8a s (P(z) + ^(5 + 6a)) 2a s v ' 



where 



z = ±\l—Z-(t-to). (3.185) 



Let us come back to the potential V Eq. (3.163) and let us look for its critical 
points. One gets 

% = — (a^ 2 - a s <P - 1) (2a s <P - 1) (3.186) 



and from 



easily finds 



dV 

- = (3.187, 



#i = (3.188a) 



<P 2 = — (3.188b) 

2a s 

<P 3 = (3.188c) 



One obtains 



25/4 

"<*.) - I^f 

y(<2> 2 ) = (3.189b) 

V($ 3 ) = 0. (3.189c) 

It is easy to see that #i is a maximum and <P 2 and ^3 are minima of the potential V. 
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However, in our simplified model of an evolution of Higgs' fields we have not a 
second (local) minimum. Thus in order to mimic a real situation we start an evolution 
of a Higgs field from the value 

3> = 0. (3.190) 

In this case 

V(0) = 4. (3.191) 

Thus 



^initial) = (3.192) 

+ \ 

2a. 



<P(t mA )= ] -^. (3.193) 



Coming back to Eq. (3.172) we get 



1 + V5 

f dip lAe n ^ . 

\ n 2 iW 2 = = V^^ ^nd -initial • 3.194 

J \/(a - if 2 + if + l)(a + v? 2 - if - 1) V r 2 
Thus the amount of inflation obtained in our evolution is equal to 

1 + V5 



^WS" / „ = , ,X , „ (3.195) 



o 



dtp 

y/(a -<p 2 + cp + l)(a + </? 2 - ip - 1) 



Supposing simply 

a < | (3.196) 

one gets 



where 



or 



^ ^ \l 71^7 / ; ^ (3-197) 

5Aae I v / h¥)^ 



cos 2 ipi = - — (3.198) 

^ 5 + 4a v ; 

cos 2 (p 2 = -^— (3.199) 
5 + 4a 



/ ^ 

JV„ = hJ r J „ ! = (3.200) 
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or 

In order to get 60-fold inflation 

N ~ 60 (3.202) 

we should play with parameters in our theory. Let us notice that in our simple model 
the constant a.\ equals 



ot\ = m 



(see Eqs (2.8-9a). 

Let us consider a slow-roll dynamic of Higgs' field in this simplified model. From 
Eq. (3.9) one gets 



Using our simplifications from Eqs (3.152-162) and Eqs (3.143-149) one gets 

d$ _ e 71 * 1 dV 
~dt ~ ~ Ur 2 H ~d$ 

supposing a condition 



(3.205) 



4 = l ffl 9pn + r'gap (3.206) 

d<P e n ^A 



v 

or 

Changing the dependent variable into 



if = a s $ (3.208) 

one gets 

dt 



UP + b (if 2 - ip - 1) (2ip - 1) = (3.209) 



where 

e n# lA 



From Eq. (3.209) one finds 



6=^->0. (3.210) 



^ " V " 1 - e -5K*-to) (3 . 211) 



2' 
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and finally 



1 ± 



-56(t-t ) 



(3.212) 



One easily notices that for (t — t ) — > oo 

i A± V5 N 



*(*) 



a. s 



(3.213) 



It simply means that a slow-roll inflation never ends. 

If we suppose that an inflation starts at <p = 3, i.e. for 



^initial 



1 , 5 



(3.214) 



it is evident that to complete it we need an eternity. 

Summing up we get in our simplified model a finite amount of an inflation for a 
nonstandard dynamic of a Higgs field and an infinite amount of inflation for a slow-roll 
dynamic. The truth probably is in a middle. 

Probably it would be necessary to consider a full equation for a Higgs field in a 
simplified model 



or 



^.-3H ^- 3bH (ai<P 2 - a s <P - l) (2a s <P - 1) = 
(Pip 



dF " 3Ho H 



3ba s H (if 2 -if- 1) (2<p - 1) 



Changing an independent variable from t to r = H t one gets 
d 2 <P dip ba s ( 2 w 



dr 2 



Let 



3ba s 



and let us change a dependent variable ip into 

z = 2ip - 1 

z + 1 
V 5 = — ^— • 



One gets 



dr 2 



dz a o 
dr 4 



(15a) 



z = 0. 



(3.215) 
(3.215a) 

(3.215b) 
(3.216) 

(3.217) 
(3.218) 

(3.219) 
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Now let us take a special value for a = — ^ and let us change z into r: 

z = iV5r. (3.220) 

One gets 

^ + 3^-2r 3 + 2r = 0. (3 . 2 21) 
This equation can be solved in general [10] 

r( T ) = ide T sn(Cie T + C 2 , -1). (3.222) 

Thus 

z(t) = -v / 5Cie T sn(Cie r + C 2 , -1). (3.223) 

sn(u, — 1) is a Jacobi elliptic function with a modulus K 2 = — 1, C\ and C 2 are 
constants. In this way 

Mt) = — (l- V5C ie Hot sn(C ie Hnt + C 2 , -1)) ■ (3.224) 

2a s V / 

However, we are forced to put a = — j^. Let us remind that 

7 A (3.225) 



272^2- 



and we are supposing A > 0. However (see Eq. (3.164) for a definition of A), it is 
possible to consider A < only for a pleasure to play Thus we use formula (3.224) to 
consider a dynamics of a Higgs field. Let us start an inflation for t = with <P = 0. 



One finds 



(i.e. initial = 0). 

Let 



<2>(0) = 0. (3.226) 
= sn(Ci + C 2 ,-l) (3.227) 



n, y 2 

(a true minimum — a "true" vacuum). One gets 

1 



#(*end) = — f (3.228) 



_ Cie ff0tend gn (C ie ^0tend + _j) _ (3.229) 
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Let us calculate the derivative of $(t). One gets 

f =-^C.^(-(Ci«*« + ft,-l) 

+ cn (C ie Hot + C 2 , -1) • dn (de" ' + C 2 , -1 



(3.230) 



Let us suppose a slow movement of <P such that 

— (0)=0. (3.231) 

From (3.231) one gets 

sc(Ci + C 2 , -1) = - dn(Ci + C 2 , -1). (3.232) 

This gives us an equation for a sum of integration constants. Thus we can get from 
Eq. (3.232) and Eq. (3.227) integration constants C\ and C 2 and from Eq. (3.229) t en d, 
which can give us an amount of inflation 

N = H t end . (3.233) 

In some sense Eq. (3.229) is an equation for an amount of inflation 

I = -C ie N " sn (c ie Na + C 2 , -l) . (3.234) 

Using Eq. (3.232) and Eq. (3.227) one easily gets 



Ci = -y 5 ^ ~ -0.5256. (3.235) 

One gets 

Ci + C 2 = -= (F(90°\45°) - F(32°\45°)) = 0.9431 (3.236) 
V2 

where F is an elliptic integral of the first kind and 

arcsin ( J 5 ~ ^ I ~ arcsin(().o2.")(j) ~ 32°. (3.237) 



V 



10 



Let us come back to Eq. (3.234) and let us denote 

e^ = x . (3.238) 
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One gets using (3.236) and (3.235) 



C 2 ~ 1.4687 (3.239) 

and finally 

1.4866x0 - 5.0354 = F(<p/45°), (3.240) 

where 

1.9026 

cos w = . (3.241) 

x 

Prom Eq. (3.240) one finds 

2 ' 8284 5.0354 / , M (3.242) 

We come back to this equation later. 

Thus we get the following results. The Higgs field evolves from a "false" vacuum 
value (<P = 0) to the "true" vacuum value (^ = <?o) completing a second order phase 
transition. The potential of selfinteracting & field changes and a new equilibrium <Pq is 
attained. However, the field \& evolves from iFq to new value a little different from <?o 
and afterwards a radiation era starts. The change of the value of a scalar field )& is 



small in terms of a variable y (see Eq. (2.69) for a definition) only from y ^5 to ^ ( n ^ s 
a natural number n = dim//, and is at least G2, dimG2 = 14). Let us consider an 
evolution of a field \P in this epoch. From Eq. (2.59) we get 



<p + w£<p = (3.243) 



2 n 

^=^(„ + 2)( r ^ ; y| T |(J2l] (3,244) 



where 



2M V n + 2 / V /3 

!f = % + ¥> (3.245) 

and we linearize Eq. (2.59) around •P'o neglecting a term with a first derivative of (p. 
Thus a scalar field & undergoes small oscillations around the equilibrium \P . These 
oscillations cannot be too long for a field should change from 



^o = ln Jr^L) (3.246) 



(n + 2)0 



to 




(3.247) 
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In terms of a field ip from ip = to ip = \ ln(l + — ^- Thus we have 

ip = (p sin(uJot). (3.248) 
For t = 0, 99 = the time Z\t to go from ip = to p = ^ is simply equal to 



arcsinf^— ) 



< 



7T 



u 2w 

Thus the full amount of an inflation in this epoch is equal to 



(3.249) 



Am = ^ 

UJ 



arcsm 



1 



(p n 



< 



7T#i 



7rM 



pl 



In our notation 



2 ^o 2v / 6M(n + 2) 
ti = t r + At 



(3.250) 



(3.25i; 



(see Eq. (2.89)). 

After a time t\ the Universe goes to the phase of radiation domination with a strong 
interaction between a radiation and a scalar field & up to a minimal temperature where 
an ordinary (a dust) matter appears. This matter evolves afterwards independently of 
a radiation and of a scalar field . The scalar field now is evolving as a quintessence. 

In order to get some comparison let us consider an evolution of a scalar field \P 
during the second de Sitter phase in a slow-roll approximation. Using Eq. (2.59) one 
gets 

dy _ By^(y 2 -(^- 2 )) 

dt 



or 



where 



and 



/ 



V^y 2 



2/^V-(^)) 
B = 2(n + 2)\j\ 



dy = B(t - t ) 



(3.252) 
(3.253) 



n 




n + 2 





Let 



/ = 



V 1 -y 2 dy 



y^\y 2 -{^h)) 



n+2' 



\A - y 2 dy 

(y 2 - (^2)) 



(3.254) 



(3.255) 



(3.256) 
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for 



n 



<y< i 



n + 2 
n > dimG2 = 14. 



(3.257) 
(3.258) 



For the integral I one finds: 

1 



I = — arcsin(y) 



+ 



- + \h: 



^n(n + 2 



1 



n + 2 



2 n + 2 



- + \hz In 



(l + y^S) ^ + 75?) + ^W 1 ^ 




(3.259) 



In order to find an amount of an inflation let us find limits of / for y = 1 and y 
One finds 



n 

n+2 • 



lim / 

s/-i 



7T 

2 



+ 



+ 



In 2 1 + 



1+2)) 



(n + 2) 



In 



*)) 



(n + 2) 



lim / = — oo. 



(3.260) 
(3.261) 



n + 2 



Thus we see that a slow-roll approximation offers an infinite in time evolution of a 



field y from y ^2 to 1. It is similar to an evolution of a Higgs field in some slow-roll 

approximation schemes. The inflation never ends. Let us come back to the Eq. (2.110) 
and let us calculate a Hubble constant and a deceleration parameter for this model 

(t_tl)3 (3.262) 



H _ R _ ( Jh\ n+4 \ 

R \3M(3 n (n-l)J (1 



+ 



RR 



1 



3M/3»(n-l) 
n+4 



+ 



, _ , . (3.263) 

Let us compare Eq. (3.262) with the similar equation for radiation dominated 
Universe in General Relativity 

2 



H, 



GR 



3(t 
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(3.264) 



and let us calculate a speed up factor 



_a_ = _, 7r< (t-t^ 26 

H GR Mfl -(„-l)( 1 + a Jgj^_ (( _ (l)4 ) 

and let us do the same for a model Eq. (3.58) (i.e. radiation dominated with a quint- 
essence). Using Eq. (3.59) one gets 

H ^R 2Btgh(H^m p i(t-t 1 )) 
Let us notice that for a (i — ti) ~ 

A x ~ (3.267) 

and 

lim Z\i = +oo. (3.268) 

t— >oo 

Thus at early stages model Eq. (2.110) is slower than that in General Relativity. 
For large time the behaviour depends on details of the theory. In the case of model 
Eq. (3.58) we have 

lim A 2 = | (3.269) 

and 



lim A 2 = +oo. (3.270) 

t— >oo 

Thus in the first case A is monotonically going from to +oo and in the second case 
from | to infinity. According to modern ideas an expansion rate in early Universe has 
an important influence on a production of light elements, the so called primordial abun- 
dance of light elements (see Ref. [11]). If at a beginning of primordial nucleosynthesis 
the Universe expansion rate is slower than in GR, then we have 4 He underproduction. 
This can be balanced by considering a larger ratio of a number density of barions to 
number density of photons (remember the barion number is a conserved quantity) 

r, = ^ (3.271) 

n 1 

where n^, n 7 — barion and photon density numbers for a high redshift z = 10 10 . In con- 
trast to the latest if an expansion rate became faster than in GR during nucleosynthesis 
process, those bigger i] (traditionally one uses the so called) 7710 

Vl0 = W 10 r] (3.272) 
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do not result in excessive burning of deuterium because this happens in a shorter time. 
The standard model of big-bang nucleosynthesis (SBBN) demands 



3 < r/io < 5.6. 



(3.273) 



It seems in the light of observational data from cosmic microwave background (CMB) 
(BOOMERANG and MAXIMA) and from the Lyman a-forest that rj 10 = 8.8 ± 1.4, 
significantly higher than the SBBM value (for CMB) and 8.2 ±2 or 12.3 ±2 for Lyman 
a-forest. Thus the second model could help in principle to explain the data without 
modification of nuclear reaction rates due to neutrino degeneracy or introducing new 
decaying particles. As the Universe expands, cools and becomes more dilute, the 
nuclear reactions cease to create and destroy nuclei. The abundances of the light 
nuclei formed during this epoch are determined by the competition between a time 
available (an expansion rate) and a density of reactans: neutrons and protons. 

The abundances of D, 3 He, 7 Li are limited by an expansion rate and are deter- 
mined by the competition between the nuclear production/destruction rates and an 
(universal) expansion of the Universe. The SBBN theory is based on a flat radiation 
dominated Universe model in General Relativity and found in a laboratory nuclear 
reaction rates. Thus it is strongly constrained. Any significant discrepancy between 
observed and calculated value of 7710 (known as baryometry) could be very dangerous. 
Thus changing the ratio of an universal expansion relative to the model of General 
Relativity can give some margin in a primordial alchemy. 

In our theory after two inflationary phases we have in principle two radiation 
dominated phases described by Eq. (2.110) and Eq. (3.58) with speed up factors (3.265) 
and (3.266). However, the important point is to match those models. We get 



R(t d ) = R exp (H t r + ffi(ti - t r )) 1 + 



(n-1) 



(3.274) 



and 



, , n+4 7min A I i-l ~r n+4 '/mi 

| 7 | 4 / \ | 7 | 4 



B = Pr (t d )R 3 (t d ) =Pr[h + 



3 

= pminRo exp (3H t r + 3Hl(t! - t r )) [ 1 + \ ^ 



(3.275) 



(n + 1) 



where p m in is given by 

Pmin = M, 2 



(3 n+l 



(1 + % 



>2(n+l) 



^ 4 (l + r/ n 



(3.276) 
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W±{y) is given by formula (2.90) and r/ min by formula (2.97). 

The time t r = t e nd in any inflation model for an evolution of the Higgs field in the 
first de Sitter phase. The time t\ = t r + At can be calculated from Eq. (3.249) 



At 



arcsinf^— ) 



LOq 



(3.277) 



and is bounded by 



At < 



1 



2M 



M p i|7|2 V (n + 2) 
From Eq. (3.274) we can get the constant to 



n + 2 



n 



l7l 



(3.278) 



t = h + 



l7l 



Tt+4 /min 



r? 2 - - T 

'mm 



(3.279) 



where 



i?arsh 



^y¥ ex P (2A"o«r + 2ffi(ti - M) (l + (^T)^ in ) 



2V3^M pl 



(3.280) 



and B is given by Eq. (3.275). 

Let us notice that for t = td the Universe undergoes a phase transition due to a 
change in an equation of state for a scalar field W from pq, = in a first radiation 
dominated phase to p# = — pq, in a second radiation dominated phase (a quintessence 
phase). 

A matter which appears in the second phase does not play any role in an evolution 
of the Universe. For t = td we have a discontinuity in Hubble constants (parameters) 
for both phases 



H{t d ) = -= 



2 | 7 |«+4 



(td-h)* 



3M(n - l)f3 r 



1 



(n-1) 'In 



(3.281) 



in the first phase and 

H(t d ) = V3m p i 

for the second one, 

1 



^ ctgh ( 



B 



B 



H(t d ) + H(t d ) 



1 2 + 1 n 2 



2 p2(n + 2)2 + 



(3.282) 

(3.283) 
(3.284) 
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Let us notice that a speed up factor A 1 changes from 

A 1 (t 1 ) = (3.285) 

to 

and an expansion rate can be slower than in General Relativity. 

The speed up factor A 2 has a more complicated behaviour than A\. It is natural 
to expect a rapid speed up of a rate of expansion resulting in non-equilibrium nuclear 
synthesis of light elements. Thus the presented scenario offers interesting possibilities 
for an evolution of early Universe. Recently some papers have appeared on scalar- 
tensor theories of gravitation exploiting the idea of a speed up factor in primordial 
nucleosynthesis [12], [13]. 

Finally let us come back to the model (3.64), (3.67). We cannot invert analytically 
the formula (3.67). Moreover taking under consideration Eqs (3.85-88) and making 
some natural simplifications we come to the following approximative formulae for R(t). 

For x < 1.1969 (Eq. (3.86)) 

/ A 

R(t) = ? — (0.44 + 0.0857V) (3.287) 
V PQ 

where 



N = exp ( 0.374^^(t - to) + 157.93 ) . (3.288) 

V m pi J 

One can calculate a Hubble parameter and a deceleration parameter and gets: 

R= R = ^ 3T7jV 

R Mpi 44 + 8.5iV v ; 

MM 

R 2 N 

We have H > 0, and q < 0. Thus the model expands and accelerates. First of all we 
consider 

p m R3 e 

and we take boundary for R. In this way one gets 

p m = a l -p Q e a Q, i = 1,2,3,4, (3.291) 

ai = 2.015 
a 2 = 0.034 

3.292 

a 3 = 1060 V ; 

04 = 8 
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If we reconsider a contemporary gravitational constant Gn as G^e ( n + 2 )* n we would 
get interesting ratios. Thus we get 

— = 0.034 -7- 2.015 or 8 -=- 1060 (3.293) 

which for the first is an excellent agreement with recent data concerning an acceleration 
of the Universe. Solving Eq. (3.288) for N if R = 1.33^J gives us 

N = 10.47 (3.294) 



H = • 0.25 (3.295) 
M p i 

-q = 1.5. (3.296) 

Using Eq. (3.288) and Eq. (3.295) we can estimate a time of our contemporary 
epoch 

t = -- ■ 10 15 yr-t (3.297) 

where h is a dimensionless Hubble parameter, H = h- Hq (we take for a contemporary 

Hubble parameter H = 100^^), 0.7 < h < 1, which seems to be too much. 
We can also estimate a density of a quintessence 

p Q = ^— . 0.571 = h 2 • 0.08 • lO" 23 -^! = 8h 2 • 10" 29 -^ . (3.298) 
v G N m 3 cm 3 

Let us come back to the Eqs (3.177) and (3.184). In Eq. (3.184) a Higgs field is 

given by a P Weierstrass function. This function can be expressed by Jacobi elliptic 
function 

P{z) = e 3 + (ei - e 2 ) ns 2 [z{e x - e 2 )*) (3.299) 

with 

ei - e 3 

where ei, e2, e3 are roots of the polynomial 

4x 3 -52X- 53- (3.301) 
Thus we should solve a cubic equation 

- 3 " 12 (5 " 3« 2 ) (t + yJ = °- (3 - 3 ° 2) 
We want Eq. (3.302) to have all real roots. Thus we need 

p = -— (5 - 3a 2 ) < (3.303) 
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and 

D = 

27 4 



3 2 

D = V - + q - < (3.304) 



where 

1 /a 2 10 \ 

Both conditions (3.303) and (3.304) are satisfied if 

0<a< 0.3115. (3.306) 



In this case we get 



ei = 1^5 -3a 2 cos | (3.307) 



e2 = - V 5 — 3a z cos — - — 



where 



From (3.307-309) one gets 



12(5 - 3a 2 ) 3 / 2 ' 



Thus 



where 



(3.308) 



e 3 = 1 ^5 - 3a 2 cos y + 4?r (3.309) 



(9a 2 + 40) 

cos</? = — — — — , . (3.310) 



9 sin -§ 

K 2 = 3.311 

sin V 7 

ei > e 2 > e 3 (3.312) 



W 8a s (^ in s 2 (n) +^ 2 ) 2a s v 7 



, Ae* (5yle n *i - 3a 2 r 2 ) / <p + vr . 
u =\l " ~4 — \/ sm ^V- (* - *o) (3.314) 



A 1 = e 1 - e 2 = \/^ - a 2 sin - ^ - (3.315) 



^2 = e 3 + -^(5 + 6a) = ^ | - a 2 cos + ^(5 + 6a) (3.316) 
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and 

yfAe^ (9a 2 Br 2 + 40Ae n ^) 

cos if = i— 2 ^-L 3.317 

12 (5Ae n ^ - 3a 2 Br 2 ) 3/2 

/ a 2 Br 2 

°<V^r <0 - 3115 ( 3 - 318 ) 

and a is given by Eq. (3.173). 

The function ns(u, K) is an elliptic Jacobi function with a modulus K given by 
Eq. (3.311) 

ns(u,K) = —±— (3.319) 

sn(u, K ) 

sn(u, K) = sinam(u, K). (3.320) 

In order to justify our treatment of Eqs (3.130) and (3.133) we consider the full field 
equations 

Rfiv ~ 2 5 = 2Au ^ u - Ag^ (3.321) 

(where A is given by Eq. (3.136)). 

However, in this case we consider A arbitrary. Let us consider static and spherically 
symmetric metric 

ds 2 = e v dt 2 - e A dr 2 - r 2 {d9 2 + sin 2 6 dip 2 ) (3.322) 

where v = v(r) and A = A(r). 

From Eqs (3.321-322) one gets 

r 2 A = e~ x (rv' + 1) - 1 (3.323) 
r 2 A = e- x (r\' - 1) + 1 (3.324) 
A' = -Av' (3.325) 

where ' means derivation with respect to r. Summing up (3.323) and (3.324) one gets 

e" A (A' + v') = 2Ar. (3.326) 

Moreover in our case A is very small. Thus we get approximately 

^ ~ (3.327) 
dr 

and 

A' -v'. (3.328) 
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In this way we go to the solution (3.137) for a small A which is our case. In this way 
A = — v and 

B{r) = e v = e" A . (3.329) 

Let us check a consistency of our solution. In order to do this we consider Eqs 
(3.323-325) in full. One gets from Eq. (3.325) 

I=e""»e"". (3.330) 

For A is very small we should suppose that no is positive and very large (/xo — ► oo). 
From Eqs (3.323-324) one gets 

^{\ + v)=2e-^ - r-e (x - v) . (3.331) 
ar 

Supposing that 

X(r )+v(r ) = (3.332) 

for an established rg > (rg is greater than a Schwarzschild radius of the mass Mq 
from the solution (3.137)). Taking sufficiently big fi we get approximately 

-^-(A + w)~0 (3.333) 
ar 

and of course 

A ^ -v. (3.334) 
Let us come back to the Eq. (3.242) and consider it for 

^ = &7r + f-£, jfe = 0,±l,±2, ... , (3.335) 

where e is considered to be small. One gets 



2.8284(-l 



,fc 



2 ° 

5.0354 = I ; d ° = + 2kK(±), (3.336) 
i ./l-I s in 2 ^ 



where 



z 

/df) 
= 1.8541 (3.337) 

is a full elliptic integral of the first kind for a modulus equal \ . For a small e we can 
write 

sine~e (3.338) 

/ -*L= -V2 £ . (3.339) 
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Thus one gets for k = 21 

2.8284 r- 

— 5.0354 = -V2e + (4/ + 1)1.8541. (3.340) 

It is easy to notice that we can realize a 60-fold inflation in our model for 

1.9026 , , . 

x = —— 3.341 

\e\ 

can be arbitrary big for sufficiently big I. 

Let us take / = 25 • 10 24 . In this case we have for e an equation 

e 2 - ({41 + 1)1.3110 + 5.0354)e + 2 = (3.342) 

e 2 - (10 26 + 6.03 54) e + 2 = 0, (3.342a) 
e 2 - 10 26 e + 2 = 0, (3.342b) 



or 



i.e. 



and finally 

|e| ~ 10" 26 (3.343) 

x ^ 1.9026 • 10 26 (3.344) 

lnz - 60.51 (3.345) 

which gives us a 60-fold inflation. 

Moreover the Eq. (3.242) has an infinite number of solutions for < ip < ^, 

ip = kir + ip, 

- 5.0354 = f d ° + 7.5164 • / (3.346) 
cos ^ J ./l-Um 2 6 



2 

with 



1 9026 

x = - -, 1 = 0,1,2,3,... , k = 21. (3.347) 

COS 'if 



One can find roots of Eq. (3.346) for 
/ = 5 

x = 16.83 
lnx = 2.82, 

/ = 50 

x Q = 0.24 • 10 3 
lnx = 5.48 



(3.348a) 
(3.348b) 
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and lnxo for I = 5 • 10 24 . In the last case one finds using 70-digit arithmetics 

In x ~ 58.479... (3.349) 

which gives us an almost 60-fold inflation. In general an amount of inflation 

7V =ln:Eo (3.350) 

is a function of / = 0, 1, 2, . . . and probably can be connected to the Dirac large number 
hypothesis. 

Finally let us take I = 10™ where n > 10. In this case one can solve Eq. (3.342a) 
and get 

2 

1 1 [(4 -10™ + 1)1.3110 + 5.0354] ' v ' ' 

Thus for xq we find 

x Q ~ 5 • 10™ (3.352) 

and 

N (n) ~ 1.60 + 2.30n. (3.353) 
In this way, for large I, Nq is a linear function of a logarithm of I, 



iVo(10) ~ 24.6, iV (20) ~ 47.6, A^ (24) ~ 56.8, 
iV (25) ~ 59.1, iV (26) ~ 61.4 



(3.354) 



or 

AT (log 10 I) = 1.6 + In/. (3.355) 

It is easy to see that Eq. (3.355) is an excellent approximation even for / = 50. 

Let us come back to the Eq. (3.29) in order to find a power spectrum for our 
simplified model of inflation. Using Eqs (3.161), (3.230), (3.235) and (3.239) one gets 
after some algebra 



P R {K) - 2.8944 ( ^ j ni ajf ^— ) CWoG ) 



(3.357) 



where 

f(x) = x ~ 2 ^ S n(1.4687 - 0.5256x, -1) 

+ cn(1.4687 - 0.5256x, -1) dn(1.4687 - 0.5256x, -1)) ' \ 
Using some relations among elliptic functions one gets 

1 2dn 4 (u, i) 

f{x) = j 3.358 

x2 (sn(u,i)dn(«,i) + v^cn(u,i)) 
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where 

u = 2.0770 - 0.7433x. (3.359) 
Let us consider a more general situation for the Eq. (3.231), i.e. 



f(0) = /> (3.360) 



where h / 0. In this way one gets 



h = C 1 F (sn(C 1 + C 2 , -1) + cn(Ci + C 2 , -1) dn(d + C 2 , -1)). (3.361) 

a* 



Using Eq. (3.227) one gets 
C 5 



V2 



— 



^( 1+ ^) 4+4 -( 1+ ^) 2 



(l+ ^ r+4 _ (l+ ^ r 

2 d<£ 



(3.362) 



sin 2 



and 



|d | = ^ 10 = . (3.363) 



v7( 1 + lr) 4 + 4 -( 1 + *) 2 



Let us come back to the Eq. (3.229) to find an amount of inflation for C 2 and C\ 
given by Eqs (3.362-363). One gets 

arccos ( ^—n— I 

V 2C 1 e N I 

^2C 2 -K(l)= I d ° -V2C ie N °, (3.364) 
or using a natural substitution 

COS ^ = 2C^' (3 ' 365) 

^C 2 -K(l) = j £ - ^L- . (3.366) 
I A/l-^sin 2 ^ 2cos< ? 
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Taking as usual 

^ = 2Ztt + § -e, 1 = 0,1,2,... , 

one gets for small e 

1 1 



e ~ 



V2(4l + 2)K(±)-2C 2 {4l + 2)V2K(±) 
(for large I) and eventually 



„ = ln (^i>) +ln , 



Let us calculate ^ for t = t* = ln( fl ^ Q ). One gets 



+m ( c '(^) + c -- 1 ) d »( c -(^) + c -- 1 ) 



Thus we can write down a Pr(K) function. 



where 



/(s) = ^(sn(Cix + C 2 ,-l) + cn(C 1 x + C 2 ,-l)dn(C 1 x + Ca- 
using some relation among elliptic functions one finds 

where ^ 

* 2 (sd(n,i) + V2cd(n,i)nd(n,i)) 2 

and 

u = v^Cix + \/2C 2 . 
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Moreover we can reparametrize (3.374-375) in the following way 



tt - V2C lX + K{\) - dV2 - j V (3.376) 

'l - \ sin 2 Lp 



o 



where 



a s h 1 Cf-25 

H " |Ci|V 5 



(3.377) 



For large C\ one gets 



|Ci| > VE. (3.378) 



u^V2d(x-l) (3.379) 



Using (3.369), 



a s h C\ 



W2K[\)l 



(3.380) 



u= ^ (x-1). (3.381) 



e 



If we take large C\ (it means, a large h) 



h=^S° (3.382) 



y/Ei 



a.. 



and simultaneously sufficiently large I we can achieve a 60-fold inflation with a function 
Pr(K) given by the formula (3.374). Large C\ means here C\ ~ 100, large I means 
/ ~ 10 25 . 

Let us calculate the spectral index for our Pr{K) function, i.e. 



s[ > dlnK 



(3.383) 



One gets 



v/2cd(u, I) - sd 3 (n, \) 



n s {K) - 1 = -2 - 2C lX \ " k m i \u U 

sd(u, ±) + a/2 cd(u, ±) nd(u, 5) 



where 

« = >/2(Cix + C 2 ) = -5^- (dJf + C 2 i?o#o) 

K 
RqHq 
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(3.385) 



dn s (K) 

A very interesting characteristic of Pr(K) is also — ^ . One gets 

am K 

dn s {K) _ ^ v/2cd(u, i) - sd 3 (u, I) 



dTnif sd(n i) + V2cd(« i)nd(n i) 

/ r- i s 1 n2 ( 3 - 386 ) 

- 2CW sd 2 («, \) - 2CW (^ Cd ^)-sd(n,,)) 

(sd(«,i) + V2cd(«,i)nd(u,i)) 

where u and x are given by Eq. (3.385). 

Let us take for a trial C\ = —0.5256 and C2 = 1.4687. In this case one finds 



v/2cd(u, \) -sd 



A) 



n s {K)-l = -2 + 1.0512a; / 1 (3.387) 

sd(u, 5) + V2 cd(u, ^) nd(u, j) 

dn a (tf) _ i Q5i2x V2cd(u, I) - sd 3 (u, ±) 



din IT ' "~ sd(n,i) + v / 2cd(u,i)nd(n,i) 

- 0.5525x 2 sd 2 (u, ±) - 0.5525x 2 (v / 2cd(n, |) — sd 3 (n, -|)) 2 

(sd(u, |) + >/2 cd(u, |) nd(u, f )) 

u = 2.0771 - 0.7433x 

K (3.388a) 

RqHq 

The interesting point is to find n s (K) ~ 1 (a flat power spectrum). One gets 

sd(u, \) + V2 cd(u, |) nd(«, §) = Cix (sd 3 (u, ±) - V2cd(u, ±)) (3.389) 

and 

sd(u, I) + v/2cd(n, ±) nd(«, ±) / 0. (3.390) 
Using (3.389-390) one gets 

^ = -2C 1 2 x 2 sd 2 (.,i) (3.391) 

if Eqs (3.389-390) are satisfied. 

In the case of special C\ and C2 one gets 

0.5525x 2 sd 2 (u ±) (3.392) 



d\nK ' K ' 2 

where x, u are given by Eq. (3.388a). 
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Let us reparametrize the Eq. (3.389). One gets 

sd(u, \) + V2 cd(u, \) nd(u, |) = - - J? 6 * 2 (sd 3 (u, |) - cd(u, ±)) . (3.393) 

For sufficiently big u one gets (the equation has infinite number of roots) 

u = u 1 + 2lK{\), i = 0,±l,±2,... (3.394) 
where u\ satisfies the equation 

i 

sn(ui,i) = ^ fl + ^^65 + V203§7- \/V20357- 65^ - 0.65... (3.395) 
Moreover for x = D K „ > we have the condition 

v 2Ci 

One finds 

sd 2 (u!,i) = 0.53... (3.397) 

Thus 

*»• (ux + -0.53 Z^+^-^A 2 



din if v 1 V2;/ " \ y^Ci / ' (3.398) 

One gets 



I = 0,±1,±2,... 



arcsin0.65 ~ 40°. 54 (3.399) 

40°. 54 

/rlf) 

= = F(40°. 54/45° )=• 0.73... (3.400) 

Taking special value of Ci and Ci one gets 

1.81 -4.68/ (3.401) 

where I = 0, ±1,±2, ... . For x > it is easy to see that I should be nonpositive. 
Practically I should be a negative integer 

I < -3 (3.402) 
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and 



fin 

s -0.29(1.81 -4.68Z) 2 . (3.403) 



dlnK 
Thus for 

Ki = R H (1.81 - 4.68/) 

we get 

n s (Ki) ^ 1. (3.404) 

The important range of In i^T , Z\ In K is about 10. 
Thus 

n s {K) ~ 1 ± 2.9(1.81 - 4.68Z) 2 . (3.405) 

Taking Z = — 3 one gets 

-300 < n s {K) < 300 (3.406) 

and 

P R {K) ~ K 1 "^^) 

is not flat in the range considered. 

Moreover we can improve the results considering Eq. (3.398) for large C\. For 
large C\ one gets 

g ~ -1. (3.407) 

Thus we find 

_^£_ „ _i. 06 ( Ul + 2Z^(I) + C x ) 2 • (3.408) 
Taking large value of C\ in such a way that 

C 1 = -2lK{\)-u 1 +e (3.409) 

where e is a small number, 

e ~ KT n , 

one gets 

d " s ~ -1.06 ■ 10" 2n . (3.410) 



dlnK 

The last condition means that we should take 

h ^ W = _^2_ ( 73 _ 2/^(1) _ 10"") (3.411) 
V5a s v5a s 

where Z is an integer, Z ~ —100. 

din 

Thus for some special values of h we can get arbitrarily small — which means 

a In A 

we can achieve a flat power spectrum in the range considered (AlnK ~ 10), 

n s (K) = l±10- (2n - 1} , (3.412) 
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n arbitrarily big. 

Let us consider the value of K coresponding to our value of n s {K) = 1. One gets 



0.73 + 21K{\) - V2C 2 

X = pr^ . 

V2C\ 

Using our assumption on a large C\ and Eq. (3.409) one finds 



(3.413) 



or (for e is small and I quite big) 
and 



x ~ 0.293 (3.415) 



K ~ 0.3R H (3.416) 
with the range A In K ~ 10. It means that 

0.3e" 10 < — — - < 0.3 • e 10 (3.417) 
RqH 

which gives us a full cosmologically interesting region 

10" 4 < -J^— < 10 4 . (3.418) 

Finally let us consider in this section two problems. First we match our solutions of 
an evolution of the Universe, i.e. both de Sitter phases, radiation dominated Universe, 
matter dominated Universe and K-essence Universe (kinetic energy dominated quint- 
essence) . In order to simplify calculations we match the first de Sitter phase to radiation 
dominated Universe using some issues from the second de Sitter phase. Secondly we 
write down solution to the evolution of Universe with a cosmological constant (from 
the second de Sitter phase), with a radiation and with a matter (a dust). All of these 
material ingredients are evolving independently in this case, similarly as a matter (dust) 
and a quintessence for the model (3.67). 

Let us consider Eq. (2.26) and let us match it to Eq. (3.58). One gets 




sh y ^^ tendj (3.419) 

R x = R e Hot ™ d (3.420) 

where 

tend = (3-421) 
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is the time of the end of the first de Sitter phase. (iVo is the amount of an inflation.) 
Simultaneously we need a local conservation of an energy, i.e. 

B 



Pr(Rl) = 



R\ 



where 

see the formula 



p r (R l ) = Z{Hl-H 2 1 ) ] 



(3.422) 
(3.423) 



Pr = A c i(!Pb) " AcoOPl) = Aco(^l) ( - 1 



H 2 o 



After some calculations one gets 



B = 



end 



( 



(3.424) 



Now we need match a solution with a radiation dominated period to the solution 
with matter dominated period (with the same cosmological constant — quintessence 
energy density). One gets 

A B 



and finally 



where 



-I 



R 2 = 0.7916 ? 



B 



sh 



'2V3^pqM p1 



B 



(3.425) 
(3.426) 

(3.427) 



(3.428) 



PQ V Po \ 
t\ is a time to match (see (3.90)). Let us consider 

iVl = H X {h - tend) 

as an amount of an inflation during the second de Sitter phase. In this way one gets 

— £ + £ (3.429) 

and finally 



A = 2.015 • 



3/4 



ars 



h 4/3 



B 



2^3 M pl ^ N 



arsh 



'lL_ sh 3/2 ( 3.43 0) 

(3.431) 
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However, we should match also a time. In this way we get from Eqs (3.287-288) 

(3.432) 



418.48 „, JVi 
to = — =M pl + — 

y/PQ H l 



or using the value pq = A c q = 10' 



-52 1 



we get 



t = 0.447 • 10 14 yr + . 

ill 

If Rq is an initial value of "a radius" of the Universe we get 

n + 2 



(3.433) 



R = R ie - N ° = 



2v / 3M ph 


/P~Q 


p r arsh (y 


Pr ) 


2v / 3M ph 


[PQ 


p r arsh (y 


[pq\ 

Pt ) 



■ cxp 



^0 + 



4 

(n+2)/8 



-iV 



(3.434) 



The only one ingredient from the second de Sitter phase is using the formula (3.428). 
Now we should match the solution (3.287-288) with conditions (3.90) to the solution 
(3.124), (3.116), (3.122). We should match a field IF, a density of an energy, a radius 
and a time. One gets 



R 3 = R exp [ ^-^ 2 - to)) = 3.07 i 1 



M, 



pi 




#(t 2 ) = — In 



25 



2n V3|7l/ "A/ 3|M| 



M pl (t 2 -t ) = ^o 



PQ + i\ = p° ein+2) *° + ^ 



where 



<5 = 2 



5M p 2 , 




3n 2 


2^/3 


n + 1 




M p 2 i 


K 



exp 



^M pl( t 2 -to) 



to /0. 



One finds: 



'3M 1 JJHn + l)y) + - h 



25 Mpi 



3 | 7 |n+l n n 



(3.435) 

(3.436) 
(3.437) 

(3.438) 

(3.439) 
(3.440) 

(3.441) 
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and t<i is large, 



From 



i? 3 = 3.072 4 3 / — 



one obtains 



(3.444) 



1 (o( o A c0 n 2 -3(n+l) ', 



From Eqs (3.90), (3.287) and (3.288) we get 

t2 = | L + ^38 Mpi (3446) 



N! 5.38 /3M 1 , / 5(n + 2) n f3 n \ . 

^ft+^^-v^^'Hv (3 - 447) 

5.38 „ 

M p i = 5.64 • 10 11 yr. (3.448) 



Finally we find Rq- 

R = R 3 exp (-^±l(t 2 -to)] = 3.072//— exp ( -^±1 (t 2 - t ) ) . (3.449) 



In this way we match an evolution of the Universe from the very beginning up to 
K-essence dominance. We find all the constants. Let us notice that the solution (3.90), 
(3.287-288) has a quite known behaviour in the theory of ordinary nonlinear differ- 
ential equations. The solution cannot be extended after some value of the dependent 
variable (and also independent). One says the solution dies. Due to this interesting 
behaviour we obtain physical consequences. Matter and a quintessence cannot evolve 
independently. We get also a ratio between a density of matter and a quintessence 
which agrees with observational data and helps (in principle) to calculate a time of our 
contemporary epoch. 
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If we use Eqs (3.287) and (3.288) we can calculate a time of our contemporary 
epoch, i.e. a time when the ratio of a density of a matter to a quintessence density is | 
(this is this ratio measured now). One gets 

^temporary = 8.97 + TT = 9-42 • 10 10 yr + ^ . (3.450) 

\/PQ H-l Hi 

In the same way using Eqs (3.287), (3.288) and (3.289) we calculate a Hubble parameter 
for our contemporary epoch. 

lOO^f 

-^contemporary ^contemporary JyJpQ (3.451) 

where 

^-contemporary — 1.21. (3.452) 

Both values are a little too big. 

Moreover one can improve the results. However, the age of the Universe is a little 
longer, 

^age of the Universe — ^contemporary ~i~ tl — 9.42 • 10 yr + 2 — — H — — . (3.453) 

Hi Hq 

In general one gets 

32.27, 3 /^- 16.38 



H ^ Q 



M p i 100 sfpcL- 227.48 
V Pm 

and 

5.21 

-q= r= + 1 

11.76, 3 /— - 5.17 

V Pm 



where is a ratio of a quintessence energy density to a matter (dust) energy density, 
or 

32.27 ?fi&. - 16.38 

h = 0.9115 . 

14.96 \ - 0.03 

V Pm 

Finally, let us express an age of the Universe in terms of the ratio . One gets 
PQ ^ = 2.667 In ( 31.87 f/^ - 14.03^ ™* + 2 W + f ^ 



or 



M — ^ = 2.82 • 10 10 yr ■ In f 31.87^ - Um) MeL + 2 (^) + (^ 



Pm) V V Pm J PQ \H\J \Hq 

7.3 



For = 0.0034 one gets 

PQ b 



t(294.12) = 14.91 • 10 10 yr + 2 (^j + 0|) 



t(0.496) = 6.76-10 10 yr + 2^^ + ' V ' 



For Bm. = 2.015 

PQ 

+ ((\ /1Qfi^ — 7fi . ID 10 wt- _l_ 9 _ _ 

At = £(294.12) - i(0.496) = 8.15 • 10 10 yr. 

The last number is a duration time of the epoch for pq is equal to the contemporary 
value of a cosmological constant. 

Let us come back to the Eq. (3.52) and consider it in the flat case K = 0. One 

gets 

R<1R ±dt. (3.454) 



PQ . E>4 i A . o , B 

3M^ L ^ 3M^ L_r 3M 2 7 

P i P i P i 



We change R into x: 



and finally get 



where 




R= a/ — x (3.455) 



0<«=|^. (3.457) 



This is an evolution of a flat Universe with radiation, matter (dust) and a density of 
a quintessence (this is a cosmological constant A c o). The integral on LHS of Eq. (3.456) 
is an elliptic integral and can be calculated. The properties of the result of calculations 
strongly depend on the value of a. Let us notice that the integral we have calculated 
here has a = (no radiation). 

The evolution of a quintessence, a matter and a radiation is here independent. 
One gets the following results for 



/ 



xdx , 

. 4 — — - • (3-458) 
+ x + a 

In general there are two cases: 

I B 2 > (3.459) 

II B 2 < (3.460) 



74 



where 

= VWEIzMzJMH (3.461) 
2^126f _ a 

and b\ > is a solution of the cubic equation 

y 3 - ay - | = (3.462) 

y = 26? (3.463) 

Eq. (3.462) can be solved by using the Cardano formulae in two cases 

1) D > 0, 

2) £>< 0, 

D = 27 - 256 ° 3 . (3.404) 

27 -256 v ; 

In case 1) one gets 



bi = \\j \J±+\\[rT- 256a 3 + ^4- i\/27 - 256a 3 > 

< a < -4t= = 0.472470393 . . . 
4^ 



In case 2) one gets 



where 



(3.465) 



6i = *7^/cos| (3.466) 



cos^ = A^ (3.467) 

_3_ 
4?4~ 



a>— , ¥»€<0,f). (3.468) 



Condition I can be transformed into 

Ab\ < a + 1 (3.469) 

for both cases 1) and 2). In case 1) this condition has no solution. It means we have 
always B 2 < 0. In case 2) we have always (3.469). It means 

B 2 > (3.470) 

for an equation 

1 cos 2 m 3 
— = hi (3 471) 

^4 cos 4 /3p 4^ cos 2 / 3 ip ' V ' ' 
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has no solution in (0, ^) and in 2tt). 
The fact that in case 1) we have always 

46 2 > a + 1 

is caused by a nonexistence of real solutions of an equation 



(3.472) 



4(a + 1) 1/2 = ^7 4 + W27- 256a 3 + yl A - \y/27 - 256a 3 



where 



< a < 



4^ 



(3.473) 



(3.474) 



Thus for sufficiently small a we always have case I (in a limit a = 0, no radiation also). 
One can express b\ and a in terms of ip: 



,4_ 3 cos 2 (§) 
1 4 v / 4cos 2 / 3 <^ 

cos~ 2/3 y?. 



4V 7 ! 



Thus one obtains in case I 



(3.475) 
(3.476) 



1 <m ) 2 + / + y ((f^f ) 2 + c ) ( ( ^f )2 + d ) (1 + c)(1 + d) 

I = — In 

2 (f5f) 2 "l 



+ Ki (b\ , a)il ^arctg 
+ L 1 (b 1 ,a)F ( arctg 



P(6i,a) 
P(b u a) 



x-P 

x — a 
x-P 
x — a 
x- (3 



,ni,gi 



Pi(h,a) 



( y/l2b\ - a - 2b\ - y^6f 



1/2 



qi 



\ y/mi - a + 2b\ + yfAb\ - a t 
>/6|&i| 



(36? + ^126? _ a ) 



1/2 



/ 26? + V4of - a + v/T26f 



3/2 



Ki(6i,a) = + v 1 v _L 

U V v^26f - a - 26? - yiftf 

2^126? - a 

ni ~ 26? + y/4b\ - a + ^1264 - a 



(3.477) 
(3.478) 
(3.479) 

(3.480) 
(3.481) 
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L 1 (b 1 ,a)= M(6i,a)H 



2(126 4 - a)V2 (y 4& 4 _ " _ ^ 



1/2 

(3.482) 



M(6 1; a) 



2b\ + ^46 4 - a - ^126? - a) (26 2 + V 4 ^ - a + ^\2b\ - a) 



1/2 



a — /3 

^(66f + a)^/46f - o + 26 2 a + ^ 1926f - 1126? a + 206f a + a 3 
+ 5766P - 2406f a + 286fa - a 3 + 66 2 ^486? - 166fa + a 2 



/ X V2 

+ (a- 46 4 )<\/126 4 -a 



x ^(126 4 - a)^46 4 - a + 46 4 ^/l26 4 - a - 46 2 a + 246? + 46 4 ^/l26 4 - a 
- y/ 1926f - 1126f a + 206 4 a + a 3 - ^ 5766f - 2406f a + 286 4 a - a 3 ^ 



(^86? ^486? - m\a + a? + 2(66 4 - a) ^ Ab\ -a - 46 4 126 4 - a 
+ 46? + 66? a + y/l92&P - 1126? a + 206 4 a + a 3 

+ ^576b\ 2 - 2406? a + 286 4 a - a 3 ^) J . (3.483) 



/ V486?- 166 4 a + a 2 + 66 2 746 4 ~^ + 26fVl26 4 -a + a ^ Q 
V 26 2 + V46f - a + V126 4 - a 

/ V486«- 166 4 a+^ + 66 2 74T 4 ^- 26 2 7i^6 4 ~^T^ ^ 



c = 



(36? + V126 4 - a) (26? - ^461 - a - Vl26 4 - a) 
V46 4 - a + 26 2 ) (V 12 &? - a - 36 2 ) 



i = ^jp±§i^p, M >, (3.486) 
V126 4 - a - 26? - V46 4 - a 
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2^Ylb\-a [yJV2b\ - a + ^/Ab\ - a - 46?) 
46? - a + 26? + ^126? - aj (^/\2b\ - a - 36?) 
2^/126? - a 

(Vl2&4 - a - 36?) (V 12& i - a - 2 &? - V 4 &? ~ «) 



(3.487) 



x 26? + J4bf - a + a/126? - a 



(2(146? - a) ^46? -a- 46?a + 2^/l926f - 1126f a + 206?a + a 3 

+ 2(6? - a)^/l26? -a- 2y / 576&P - 2406f a + 286?a - a 3 ) . (3.488) 



Case II 



1 e(|Ef ) 2 + / + 2 J((fEf ) 2 + c) ((fEf ) 2 + d) (1 + c)(l + d) 
/ = - In — 

2 (f5f) 2 "l 



+ K 2 (bi,a)II I arccos 



+ L 2 (bi, a)F ^arccos 



P 2 {h,a 
P 2 (6i,a) 



x-/3 

x — a 
x-f3_ 
x — a 
x-13 



,n 2 ,<te 

12 



(3.489) 



Q2 



n 2 



1 (36? + V126? - a) 



1/2 



3i 



26? + ^46? - a + ^126? - a 
2^126? - a 



(3.490) 
(3.491) 



1/2 



f 26? + ^/46? - a - ^126? - a) ^126? - a f ^46? - a - 6i) 
^ 2 (6i,a) = ^ ; ^ — ^ TJJ^— (3-492) 



(Vl26? - a + 26? + ^46? - a) (^126? - a - 36?) 



1/2 



2V ^ ^126? - a + 26? + ^461^/ 
L 2 (6i,o) = M(6i,a) 

1 /2 

4(126? — a) (y46? - a - 6?) (V 12b i - a - 26? - V46? - a) 



(36? - ^126? - a) (26? + ^/46? - a + ^126? - a) 



(3.493) 



(3.494) 
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The function (3.456) cannot be inverted globally. It can be inverted only locally 
in some intervals where the solution lives. The solution dies on the end of an interval 
being reborning on beginning of a next interval. All the intervals can be obtained from 
the condition 

W > (3.495) 

where 



W = 



e(fEf ) 2 + / + V((f5f ) 2 + c) + d)(l + c)(l + d) 

I x-a \2 _ 1 
\x-B> 1 



(3.496) 



and 



x — a 



x-P 



/I 



(3.497) 



In the case I c > 0, d > 0. 
In the case II c > 0, d < 0. 

Thus one can get very interesting behaviour from the physical point of view, 
because every end of an interval of living (existence) of the solution of an evolution 
equation means a start of nontrivial interaction between a matter, a radiation and a 
quintessence. In all of these intervals one can find an approximate inverse function, i.e. 
one can write 

R = F(t) (3.498) 



where 



and 



F(ti) =R 1 <R<R 2 = F{t 2 ) 



h<t< t 2 - 



(3.499) 
(3.500) 



Moreover we do not develop this project here in details. F and 77 are usual elliptic 
integrals of the first and of the third kind given by the formulae (3.71) and (3.73). 
For future convenience of this project we find roots of the polynomial 



P(x) 



x 



+ x + a, a > 0. 



(3.501) 



First of all according to the Ferrari method we should solve a resolvent equation which 
is a cubic equation. One gets 

z 3 - 4az - 1 = 0. (3.502) 
The discriminant of Eq. (3.502) reads 

27 - 256a 3 



D 



108 
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and we have two cases 



D > 0, < a < jy=, (3.503) 
D<0, a>^=, (3.504) 



In the case (3.503) one gets 



3/^/27 - 25603 + 3^/3 3 / v / 2T^25fc 3 - 3^/3 , Q cnc , 

Z1 = V 67! V 67! (3 - 505) 



z 2,3 = 7: 



1 / 3 / V27 - 256a 3 - 3^ 3 / v / 27^256a 3 " + 3^3 

2 I V 67! V e7! 



z 3 =4 A /^cosf^ + ^ ) (3.509) 



± iV3 I 3 / \/27 - 256a 3 - 3^ | 3 / V27 - 256a 3 + 3^ \ 5Q6) 

2 I y 6\/3 v 6^/3 ; 

In the second case (3.504) 

Zl =4^| cos (|) (3.507) 

--VMM) (3 - 508 » 

a (Tp 47T 
, 3 C ° S l3 + T 

where cos 7p = ^ . 

We can also distinguish the special case D = (a = j^j): 

= \/I (3.510) 

z 2)3 = -±#4 (3.511) 

In the first case we get one real and two complex conjugated roots, in the second 
three real roots, in the third case two different real roots (one of them is double). 
According to the Ferrari method one gets 

2x! = ^T + ^2 + ^ (3.512a) 

2x 2 = y/z{ - ^fz2 - y/z3 (3.512b) 

2x 3 = -^+^-^ (3.512c) 

2x4 = -^+^+^ (3.512d) 
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ZlZ 2 Z 3 = 1. 



(3.513) 



In the future analysis of an evolution of our model of the Universe there are three 
important cases for z\, z 2 and z 3 : 

A. All zi,Z2, Z3 are real and positive. In this case all x\, x 2 ,x 3 and X4 are real. 

B. One root (e.g. z\) is positive, two remaining (z 2 ,z 3 ) are negative. In this case we 
have two pairs of conjugate roots. 

C. One root (e.g. Z\) is positive, two remaining (22,23) are complex conjugate. In 
this case we have two real roots and one pair of complex conjugate roots. 

This analysis is quite important because only for positive real roots the integral / 
can have singular points. It means in A and C case. This gives us some restrictions on 
the parameter a. 

Let us notice that in the case of a = (considered before) we have to do with the 
case C. Fortunately one real root is zero and the second real root is negative (x 2 = — 1). 

Let us notice that our case with a = is in some sense exceptional from the point 
of view of a general theory. In that case some of the formulae are singular and because 
of this it was considered separately. 

Let us notice that for sufficiently big a the polynomial (501) has no real roots. 
This is true for a > -^=. In this case the integral (458) has no singular points. This 
is case II. 

Conclusions 

In the paper we consider some cosmological consequences of the Nonsymmetric 
Kaluza-Klein (Jordan-Thiry) Theory. Especially we use the scalar field \P appearing 
there in order to get cosmological models with a quintessence and phase transitions. 
We consider a dynamics of Higgs' fields with various approximations and models of 
inflation [14]. 

Eventually we develop a toy model of this dynamics to obtain an amount of infla- 
tion and Pr(K) function (a spectral function for fluctuations). We calculate a spectral 
dn 

index n s (K) and — — s — for this model. We calculate a Hubble parameter and an age 
a In K 

of the Universe. 
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